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Section 1 


Two Dimensional Subsonic Flow 


of Compressible Fluids 


Summary 


The basic concept of the present peper is 
A 
to use tee tangent line to the adiabatic pressure- 


-volume curve as an approximation to the curve 
itself. First, the general characteristics of 


such rg fluid are shown. Then in Section Ia 


asf td = 


theory is developed Which in is 
l similar to thet of DemtchKÉfenko and Busemann 
but is more general and can be applied to flow i 


with velocity approaches! thet of sound, |] The H ! 
theory is then applied to calerlate the flow KÉ 
; over elliptic cylinders. In Section II the n. 


work of F. Bateman is applied to this approximated 


l: adiabatic fluid and twa equations are given which OG n ^ 
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TWO-DIMENSIONAL SUBSONIC FLOW 


OF COMPRESSIBLE FLUIDS 
by 


Hsue-shen Tsien 
California Institute of Technology 


Introduction 


Assuming that the pressure is a single-valued function of 
density only, the equations of two-dimensional irrotational motion 
of compressible fluids can be reduced to a single non-linear equation 
of the velocity potential. In the supersonic case, that is, in 
the case when See flow velocity is everywhere greater than that of 
local sound velocity, the problem is solved by Meyer & Prandtl and 
Busemann using the method of characteristics. The essential difficulty 
of this problem lies in the subsonic case, that is, in the case when 
the flow velocity is everywhere smaller than but near te the local 
sound velocity, because then the method of characteristics cannot 
be used. Glauert & Prandtl (Ref. 1) treated the case when the 
disturbance to the parallel rectilinear flow due to presence ofa 
solid body is small. They were then able to linearize the differential 


equation for the velocity potential and obtained an equation very 


similar to theme for the incompressible fluids. But there are usually 


stagnation points either in the surface of the body or in the field 
of flow, where the disturbance is no longer small. Hence, it is 
doubtful whether the linear theory can be applied to the flow near & 
stagnation point. On the same ground, the theory breaks dom in the 
case of bodies whose dimension across the stream is not small compared 
with the dimension parallel to the stream, 

To dS Bruins isch, Janzen 


and Rayleigh developed the method of successive approximations. This 
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wh 
method ma—ammis-neid0pizzskentiz-ami/put into a more convenient form 
by L. Poggi (R&f——2j and P.A, Walther -(Re®=3}, Recently C. Kaplan 


(Ref. 8) treated the case of flow over Joukowsky airfoils and elliptic 


cylinders; using Poggi's method, However, the method is rather tedious 
and the convergent very slow if the Seet, velocity of sound is approached. 
Pyar Molenbroek (RS) and Tschapligin (REf—5) suggested ^ 
the use of the magnitude of velocity # and inclination £ of j A 
A chosen axis as independent variables, and were thus able 
to reduce the equation of velocity potential to a linear equation. 
This equation was solved by Tschapligin, (ESZv—9) emicrecentiy—put. l 
into-e-mere—cenverrterrt form br-Fr—0Sieuser—and-ilr—Gienmer Lët, H 
The solution is essentially a series sach term of which is a product = 
of a hypergeometric funotion of W anda trigonometric function 
of Bf . The main difficulty in practical application of this 
solution is to obtain a proper set of boundary conditions in the 
plane of independent variables wp ang put the solution in a 
closed form. 

Tschapligin a in that a great simplification of p 
the equation in the hodograph plane results if the ratio of the specific 1 
heats of the gas is equal to -1 . Since all real gases have their 
ratio of specific heats between 1 and 2, the value -1 seems without h 
practioal significance. It was Demt@chenko (Ref. à and Busemann a 
(Ref. 4) who Log lear the meaning of this specific value of A. d 
They found that this really means to take the tengent of, pressurs- 
volume curve as an epproxinaté to the curve itself. However, they apt ta d 
limit themselves to use the tengetat the state of semt—of the ges, wu de eg 


Thus their theory oan only be applied to flow with velocities up to gel jet 


dy 
one half TUUS US Ra deed Recently, during & personal 
discussion, Th. von Kármán suggested to the author that the theory can 
be generalized &o usa the tangent at the state bs corresponding 
to the undisturbed parallel flow. Thus the range of usefulness of 


sei theory can be greatly extended, This is carried out in the first 
Melita 


past of the present EE TA 


NENNE theory, based upon 
Demtpchenko and Busemann's work,is applied to the case of flow over 
elliptic cylinders and the results compared with those of S,G. Hooker 
(Ref. 39) and C. Kaplan (Ref. $. Furthermore, results calculated by 


Glauert!w-Prandtl\s linear theory are also included for comparison. 


e 
Recently, H. Bateman (Ref. 4) demonstrated a remarkable 


EE 
reciprocity ef two fields of flow of two fluids related by a certain 
Seema Sechon) 
point transformation, It will be shown in the thied-pert of this 
[o 
paper that the flow of incompressible fluid and the flow of compressible 
the 
fluid approximated by the use of tangent to adiabatic pressure-volume 
curve can be interpreted as such a point transformation. It is thus 


possible to obtain a solution for compressible flow whenever a solution 
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are known, The characteristic’ of the 


incompressiblg flow can either be obtained by the well-known method of 
otical 
teristics at 
ly 
lieved to 


[ Approximabop +o the Adiabahe Relation | $ | 
ER i S E 2 H 


RE ar 
H , If f is the pressure, // is the specific volume and bal 


€ is the ratio of specific heats of a gas, the adiabatic relation 


T 


b TAL = constant is a curve in the pv plane as shown in 


d i) 
b Fig. la, Taw conditions near the point / f, U, ) which corresponds 
d 
! to tie state of Fruit undisturbed flow can be approximated by the 
tangent to the curve at that point. The equation of the tangent at 4 


= 


this point can be written ası 


hf = OC: else el 
ER aes um te E ftmt (a mi 
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(1) 


where 1% is the slope of the tangent and B 4 is the density of the 
. won The slope Ż must be equal to the slope of the curve at 
the point ( KA ) » therefore, 


(= (#)- GZ. IT ae 
where 4, is the BER ae 1 the conditions j, L- 


Thus eq. (1) can be written as‘ 


j peque) 


spuni makum bu 
This is an appreximate pressure-density\toy&diabetic) relation, and 


(2) 


is shom in Fig. Ko rue adiabatic relation. 


The€*Bernoullig^ theorem for compressible fluids is 
Jj 
fa- pa- f 4 (5) 
2 


where “J is the velocity of the gas and the subscripts 2 and 3 
denote two different states of the fluid. By substituting eq. (2) 
into eq. (3), the following relation iscbtained: 


x d 
ta' ta'- ies] s Si 


Nowif 4,70 ung ma Sass 5 with 
the subscript Ô denote the state ot) mest, eq. (4) gives: 


(4) 


If the square of sound velooity 4 is defined (as 4d 
usually done) as the derivative of Á with respect to f x 


eq. (2) gives: 


£y ei 2 
a = fir- 2^ = constant (6) 


Therefore, eq. (5) can be written as: 
(7) 


(8) 


It is interesting to a that from eq. (8) the density 
decreases BEE increases, ee Thus eq. (6) eee at 
the local velocity of sound increases as the velocity increases, This 
is just opposite to the real gas, because in the oase of an adiabatic 
flow of a real gas it is well known that the temperature ot ea 
decreases as the velocity Gili re enè thus the local 
sound velocity]also decreases, However, in the present approximate 
theory, the ratio z or Mach's number still increases as the velocity 
increases, as can be seen by eq. (7). But this ratio only reaches the 
value unity when e —0 , or from eq. (8) wen sv = e, It is 
thus seen that the entire regime of flow is subsonic and thus the 
differential equation See potentiel/e£-an-innempunsathie. This 
is the reason why the complex representation of, velocity potential and 
ee funetion is possible for all cases, as will be shom in he 
following paragraphs. However, one should realize that the portion of 4+ 
tengent that could be used as an approxi mation to the true adiabatic 


relation is that portion which lies in the first quadrant. Thus the 
upper limitler velocity,for practical application of the theory e OCC YQ 


po . By using eqs. (17) and (18), this upper limit 
is found to be 


(hue By l-lif, 


Since the point ( bs) lies on the trwe-edtevetic curve, the 
relation 4, = #1 can be used, and eq. (9) becomes: 
^ 


M i 
= f EP FH dE 37 An = 
This relation is plotted in Fig. 2ff Since for most practical cases ee Za See 
it is not likely that the ratio bi will rise to values much 

higher than 2 , f will remain positive, and this theory Mn give 

an approximate solution, , 


If the flow is irrotational, there exists a velocity 
potential d such that 


Q1) 


1 


where U, V are the components of 4/ in, "y and ER 
respectively. The equation of oontinuity, 


A (EO $9) = 


will be satisfied, if the stream function 4l is introduced such 


Ze? -$:-£É 


that 


Now if the angle of inclination of the velocity 4 to 
the Y axis is 8 ; eqs. (11) and (12) give: 


df - ua Ay + wt gra? 

df=- wier dal ph (13) 
Solving for fy and CS? | 

d: Sf dp EEG AN 

y- pdt we 


So long as the correspondence between the physical -ghume and hodograph 


(14) 


plane is one to one, or mathematically DE #0 » Ona—cun 
Ru hu expressed Y and ji as functions of Ae oy, and y 


as functions of W, . Thus, 
dp d du * $4 (15) 
A = Var dw + V, s "n 


where primes indicate the derivative with respect to variables U 


indicated as subscripts. Now substitutiweg. (15) into eq. (14), cna has : 7 fu t 7 


b= (fg, E, kai (PEK FEN 
H(t ES u iuo 


Since the left-hand side of eqs. (16) are exact differentials, eme—een 


apply the reciprocity relation T. akiert 


An S o) al N) 
en). lan) " 


the aid of eq. (7), eq. (17) gives: 


nr Yen Ae EES 


Zi, 5 (18) 
Ed" re, 
Solvin So x 
g for = and ds 5 
A m EP Ed 5 
"b ; wo A" 
i * E £ E ER (A (19) 
EE. Eq. m can A further simplified by introducing a 
new variable CO , such that 
B. E dae. 
Then eq. (19) becomes: / H 
GER 
" 7] (21) 
| = Vo 


| This can be easily recognized as de RiemanntCauchy Na differential 
equationw)and thus $ 44 y must be an analytic function of 


w- i . However, for convenience of calculation, another new 


set of independent variables l= W eo 4 , 
y- W tev are introduced where W= ae” 


Then eq. (21) can be written as: 


Ae, E ata 
| "E ou AC V) m 
pt 24. . af 
| atv) ou 
^" aan by integrating eq. (20), 
Ag d W- -= er (23) 


4, x +4, 


Carrying out these differentiations and simplifying with 


t 


dau (24) 


EE Is 


: € 208 
Substituting into eq. (8), were dy eme par few de desks, abo e € 
g reek 
(25) 
Eqs.(22), (23), (24) and (25) are the basic equations of 


the present theory. One recognizes eqf. (22) as the RiemanntCauchy 
differential equationġ, and thus the complex potential LE = g fi 


must be an analytic function of Die = -iV por: 

guf- F/M-W)- Fly) 
p-ip- FAV) - Ei 

Now, it is necessary to find the values of L nae d 


om - i d^. 
es st A a given set of values of M and V 3 1.9.5 Ni l 


i to find the transformation from hodogreph plane to physical plane. 
i 


By using eqs. (24) and (25), eq. (14) can be written as: 


l A: Ai As UM rds (27) D 
li ge aol tr from m 


L $ 24:2 
where W = A + H . These equations can be combined into 


| m 
E one equation by means of eq. (26). Thus, d 
T = a. 
N AE 7 
H dp = Au +i xta = EE. (28) A 
$ W ^, 
| ur — 
Hence, if an analvtic funetion VA fi Ww) is given for each value of RK 


d W , the corresponding real velocity Af can be calculated by 


eq. (24). Then the coordinate of the point in the physical plane at 


which this velocity occurs can be calculated by means of integrating 
eq. (28). The pressure at this point is then given by eq. (2). However, 
there is an objection to this procedure. That is, the investigator 


au Zei f) UT, e éi 


z Perte whether wish the chosen function //1/7), we 
Fp MS ae dS the desired shape ot" żolid boundary and flow pattern , 
Veet. In other words, this procedure still suffers 
the difficulty of boundary conditions as D is common to all hodograph 
methods, Transfermabon from Iwompresible How to Compressible Haw B men v 
However, due to the particular simple relation of eq. (28), 


ome ascertain approximatel e result shape of the by 


starting with the psettemter function, 


FÍW)- did - W Gle) 


(29) 


where W, is the transformed undisturbed velocity to be interpreted 


by eq. (23), and Ç is the complex coordinate p . This 


function is so chosen as to give the flow of incompressible fluid over 


the desired body shape in coordinates y and 7 > i x | 
dw. E De ELS 

Fre real velocity u plene of the incompressible fluidas 
v 

transformed velocitv W in the hodograph plane RT EE 1 


fluids Wen It is known that 


(30) 


(31) 


SCC. EN eqs. (30) and (31), inte eq. (28), Kammer dai E 
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A, RO: 
SS um - 
"d 4 re f dc 2/4 


Integrating; 


4 ar A 
er (E = 


Therefore, the complex coordinate in the physical plane 
Min fluid is equal to the corresponding complex coordinate 
in the physical plane ofYincompressible fluid plus a correction term. 
Since this correction term is usually small, the resulting shape of the 
body will be quite similar to the one infincompressible fluid. The 
factorjin the correction term depends upon he Mach's number of the 
undisturbed flow only. This can be saan by means of eqs. (7), (8), 


and (23), because from those equations ems—hme the following relations? m PR 


44; 2 
dw ZA Ce o OS 
« E e EZ (33) 


ANZ 
where GJ is the Yach's number of the undisturbed flow. A ale. 
8 
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first obtained from eq. (30) and then with eq. (23): 


Ying Cé 
[A apr re (34) 
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In-thia-TERb, the theory dexekeped ims—Part—l will’ be 


applied to calculate the flow over an elliptic cylinder at zero angle 
of attack. The incompressible flow over an elliptic cylinder in the 
complex coordinate C, can be obtained by applying Bee Joukowsky's 
transformation to the flow Mtn eircular cylinder in the complex 


(AA 
coordinate E with Ee at the ofiginv Therefore, 


the function FI) or W,6(2) can be written as: 
Lt [e+] 


(36) 
= 2 
Few, [e+ 4 

m nA on qud (87) 


A 
It is canvenient to carry out the computation by using the G 


coordinates. Thus eq. (32) is rewritten in the following form: 
de 
(er 2) HE) fe (38) 
Géi 
C A 
If conditions over 


the surface of the elliptic cylinder emiy, Auer emeerrud, Au. 


Zeta" 


ġew 


(39) 


where b is the argument as shown in Fig, & and b is the radius of 

h , 

i the circular section in the Č — plane which determines the thickness 
ratio of the elliptic section in the € - plene. Substituting eqs. (36), 


(37) and (39) into eq. (38), and carrying out the integration, the 


H following expressions for the ¥ and ri coordinates corresponding to 

Iv & d are obtained by separating the real and imaginary parts, 

I". Sax A? Ez SW ZA 
j 57 ge(pad)ont Y {Mie ewe GY pg HN 
, V (1*1) Ja y 4 KE (434 cob) 


y D [1- 7) T Ay i 4 (I-E) tap4 (21,7 A E 


The horizontal and vertical semi-axis of the elliptic section can 


then be calculated by substituting f=9 and S=F 


respectively into eq. (40). The thickness ratio d is thus obtained 


| el t HER EE 
| (BA T EE © 


For a given thickness ratio and Mach's number ER E 


the value of TM) is first computed by means of eq. (33), and 
fi k 
then eq. (41) is solved graphically for D 2 
| 
i: After "1 is obtained, the coordinate E and P for 
| 


| ' each value of H can be computed by using eq. (40). It is fortunate 
that the values of X, / so obtained lies very close to the true 


| elliptic section nn Me Hence, the velocity and the 


pressure distribution obtained by using eqs. (34) and (35) are considered 
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Critical Velocity Jor Fh; he l- Iinders 
If the velocity of flow over a body is gradually increased, 


the maximum local velocity in the field will also be increased. When 


uf 
a : the maximum 


local velocity reaches the local velocity of ae —Them shock waves 
RR on drag of the body suddenly increases,  Themméeme, 
eis velocity i ar sone. derable interest to practical engineers and 
is usually referred to as the critical velocity of the body. It is 
shown by Kaplan (Ref. $) and others that at this critical condition 


the ratio of maximum velocitv M may in the field tofthe undisturbed 


TERE YAT ern 


velocity W, , is related to WKÉ Mach's number AL bc P 


2 


fru Simemm in the following manner: 


£ 
d, Je ër: s uw 


Way, in the flow over an elliptic cylinder at zero angle of 


attack occurs at the top of the cylinder. Using eqs. (34) and (33) 
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Equating eqs. (42) and (43), the equation for calculating 
the critical Mach's number (4b Jet the undisturbed for 


each value of b is: 
sch 
it my o EE 
Ate" Meet ri 


gi Knowing P for each value of b » the 


corresponding value of fu can be calculated by means of eqs. (34) 


ae end (41). re A shows the result of this calculation with Kaplan's d 


(a 3 er 


palin PD, 
ar d Date 


a eil eu 


value inoluded for We dom It is seen that the critical Mach's 


inte, a mre firtina) 
number is lower than these obtained by Kaplan. This er 


effect of compressibility auta &hom--imbiucummendéeéhenry and is 
M 
! consistent with the results shown in Figs, ri + H n 
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If two new functions x and Y are defined by 


fal = fly tut 


JU 3 (45) 
au RAO LO csi uer duet inal, a 
drag-and Lift £funebtons. Bug means of eqs. (11) and (12), it oan 
be shown that the following relations hold: 
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(49) 


"EE, The quantities X ana SN have the dimension of a 
velocity and nn ag can be considered SAME 
in the plane whose coordinates are denoted by L and E . This 
relation between the L rá plane anà the Z ard y plene is shown 
in Fig. Ay. is thus seen that if the O REMTE E the 


SC? plane is in the positive X -direetion, the undisturbed stasom A / 


in the FE plane will be in the negative T direction. Furthermore, 


I Ben there exists a complete reciprocity between the Ra plane and 
i e 
1 the / Y plane, as shown by H. Bateman (Ref, 4). 


` 


< — ——- Comparing eq. (51) with eq. (12), it is evident that ĠU 
can be considered as the density of a fluid in the new X Y plane. 
maken sf with Incompreikje Flaw 
o far the relations obtained are general, i.e., they 


apply to fluids of arbitrary properties. However, since only the flow 
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M 
of incompressible fluid » it would 5e 
M Aa to DEEP LV olme it 
the shew in the Æ 4 plane is incompressible. If the fluid in the 
L 4 plane is incompressible, then = = / , end the Bernoullitg — Z 
hea : Heorem ges E 


S 
+ Erde Sé 


(52) 
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denote the pressure)and e denote 
2 3 
VA % A 3 th rnoulli%g theorem gives: 


PE ET constant (53) 


«—— — In view of eqs. (49), (50), and (52), eq. (53) can be 
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written in the following form: 


$ et PEE p KJA constant (54) 


-————————— By differentiating eq. (54) with respect to Z a 
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multiplying the resulting expression by "o and then integrating 
o 


~ 
with respect to E , Gas-àes the following relation conneoting 


the pressure P and'densitv 0” for the fluid in the T plane xo ef 3 


To ei iat C -4r il (55) 
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A -————— Comparing eq. (55) with the approximate adiabatic relation 


eq. (2), also noting eq. LE it is evident that eqs. (55) and (2) 


are identical If. 


In eq. (56) A is the velocity of sound in the P2 46577 and the 


NA) l refers to the conditions in the undisturbed pamm 
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TWO DIMENSIONAL SUBSONIC FLOW 


OF COMPRESSIBLE FLUIDS 
by 


Hsue-shen Tsien 
California Institute of Technology 


Introduction 


Assuming that the pressure is a single-valued function of 
density only, the equations of two-dimensional irrotational motion 
of compressible fluids can be reduced to a single non-linear equation 
of the velocity potential. In the supersonic case, that is, in 
the case when the flow velocity is everywhere greater than that of 
local sound velocity, the problem is solved by Meyer ^ Prandtl and 
Busemann using the method of characteristics. The essential difficulty 
of this problem lies in the subsonic case, that is, in the case when 
the flow velocity is everywhere smaller than but near to the local 
sound velocity, because then the method of characteristics cannot 
be used. Glauert & Prandtl (Ref. 1) treated the coase when the 
disturbance to the parallel rectilinear flow due to presence of a 
solid body is small. They were then able to Jinearize the differential 
equation for the velocity potential and obtained an equation very 
similar to those for the incompressible fluids. But there are usually 
stagnation points either in the surface of the body or in the field 
of flow, where the disturbance is no longer small. Hence, it is 
doubtful whether the linear theory can be applied to the flow near a 
stagnation point, On the same ground, the theory breaks down in 
case of bodies whose dimension across the stream is not small compared 
with the dimension parallel to the stream, 

To treat cases in which the body is blunt nosed, Janzen 


and Rayleigh developed the method of successive approximations. This 


method was explained physically and put into a more convenient form 

by L. Poggi (Ref. 2) and P.A, Walther (Ref, 3). Recently C. Kaplan 

(Ref. 4) treated the case of flow over Joukowsky airfoils and elliptic 
cylinders, using Poggi's method. However, the method is rather tedious 
and the convergent very slow if the local velocity of sound is approached. 

Molenbroek (Ref. 5) and Techapligin (Ref. 6) suggested 
the use of the magnitude of velocity W and inclination Ø of 
velocity a chosen axis as independent variables, and were thus. able 
to reduce the equation of velocity potential to a linear equation. 

This equation was solved by Tschapligin (Ref. 6) and recently put 
into a more convenient form by F. Clauser and M, Clauser (Ref. 7). 
The solution is essentially a series each term of which is a product 
of a hypergeometric function of # anda trigonometric function 
of f . The main difficulty in practical application of this 
solution is to obtain a proper set of boundary conditions in the 
plane of independent variables A. / and put the solution in a 
closed form. 

Tsohapligin (Ref, 6) shows that a great simplifioation of 
the equation in the hodograph plane results if the ratio of the speoifio 
heats of the gas is equal to -1 . Since all real gases have their 
ratio of specific heats between 1 and 2, the value -1 seems without 
practical significance. It was Demtschenko (Ref. 8) and Busemann 
(Ref. 9) who made clear the meaning of this specific value of -1, 

They found that this really means to take the tangent of pressure- 
volume curve as an approximate to the curve itself. However, they 
limit themselves to use the tangent at the state of rest of the gas. 


Thus their theory oan only be applied to flow with velocities up to 
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one half that of sound velocity, Recently, during a personal 
discussion, Th, von Karman suggested to the author that the theory oan 
be generalized to use the tangent at the state of gas corresponding 

to the undisturbed parallel flow. Thus the range of usefulness of 
the theory oan be greatly extended, This is carried out in the first 
part of the present paper. 

In the second part of this paper this theory based upon 
Demtschenko and Busemann's work is applied to the oase of flow over 
elliptic cylinders and the results compared with those of S.G. Hooker 
(Ref.10) and C. Kaplan (Ref. 4). Furthermore, results calculated by 
Gleuert's-Prendtl/$ linear theory are also inoluded for comparison. 

Recently, H, Bateman (Ref. 11) demonstrated a remarkable 
reciprocity of two fields of flow of two fluids related by a certain 
point transformation, It will be shown in the third part of this 
paper that the flow of incompressible fluid and the flow of compressible. 
fluid approximated by the use of tangent to adiabatio pressure-volume 

curve oan be interpreted as such a point transformation. It is thus 
possible to obtain a solution for oompressible flow whenever a solution 
of incompressible flow is known. The difficulty, however, lies in 

the proper choice of this incompressible flow so that the solution of 
compressible flow will have a solid boundary desired. Fortunately, this 
transformed boundary is independent of the Mach's number of undisturbed 
flow. Therefore, by a modification of the procedure, two relations 
between the velocity and the pressure distributions over a body at low 
velocities and those at high velocities over the same body are obtained, 
This enables one to predict the characteristics of the flow over a 

body at any Mach's number whenever the low speed characteristics of 


the flow over the same body are known. The characteristics of the 


incompressible flow can either be obtained by the well-known method of 
conformal mapping or by experiments. Due to the fact that practical 
aerodynamic engineers usually have the low speed characteristics at 
hand and that high speed data have to be obtained by use of a costly 
high speed wind tunnel the above mentioned relations are believed to 


be of considerable use to them, 


In the fourth part of the paper, the theory developed in e 
Part III is applied to oorrelate airfoil data obtained by J. Stack (Ref. 12) l 
in the N.A.C.A. 24" high speed wind tunnel, The agreement with theory 
is found to be satisfactory. Then this theory is applied to predict 
the compressibility effect on the lift and moment of N.A.C.A. 4412 
- A. airfoil using experimentally determined pressure distribution over 
the same airfoil at low speed. The result is again compared with the | 


more simple Glauert-Prandtl theory. 


| f | 


If f is the pressure, V is the speeifio volume and 


7 is the ratio of specifio heats of a gas, the adiabatic relation 
£y! - constant is a curve in the fr plane as shom in 
Fir. la. Now oonditions near the point f 4, V J which corresponds 
to the state of fluid at undisturbed flow can be approximated by the 
tangent to the curve at that point. The equation of the tangent at 
* this point can be written as 


p-p- Giele Ge e) 


(1) 


] 
4 
d 


wiere E is the slope of the tangent and € is the density of the 
gas. Now the slope e must be equal to the slope of the curve at ch 
the point f, y, » therefore, i 


d E) E SE Je mn. 
where (2 is the sound velocity kd żero, À the conditions 4.1. d u ; 


Thus eq. (1) oan be written as 


DE PNE DI es ES m. Li 
Á f 4, S, ($ £/ (2) i 
This is an approximate pressure-density to adiabatio relation, and Kä 


is shown in Fig. 1b with true adiabatic relation. 


The Bernoullis' theorem for compressible fluids is 


iw-fa- /4 (3) 
A 


where ww is the velocity of the gas and the subscripts 2 and 3 
denote two different states of the fluid. By substituting eq. (2) 
into eq. (3), the following relation isdtained: 


ET AI ee n 
Coo ciet ua Mean. eel 5. (4) ven 


Now if Ww =) F WW, $ ven P , and i =f » with ty 
the subscript CO denotes the state of rest, eq. (4) gives: e 


am d 
WES aUe. ute 


bAa ki 2 (5) 
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EL 
If the square of sound velocity @ is defined (as 


usually done) as the derivative of D with respect to Ba ` 


PT eq. (2) gives: 


i^ fi bra 
AR tel- Be. 4,8 (6) 


de f = constant 


Therefore, eq. (5) can be written as: 


9 2 

($J -/- 4 Ge 

Similarly, fos aw (8) 
i ee 


It is interesting to notice that from eq. (8) the density 
decreases as velocity increases, as expected. Thus eq. (6) shows that 
the local velocity of sound increases as the velocity increases, This 
is just opposite to the real gas, because in the case of an adiabatic 
flow of a real gas it is well known that the temperature of gas 
decreases as the velocity of gas is increased, and thus the local 
sound velocity also decreases, However, in the present approximate 
theory, the ratio zr or Mach's number still inoreases as the velocity 
increases,as can be seen by eq. (7). But this ratio only reaches the 
value unity when © =0 „or from eq. (8) when éi en, It is 
thus seen that the entire regime of flow is subsonic and thus the 
differential equation of velocity potential/efen—incompresstbie., This 
is the reason why the complex representation of velocity potential and 
stream function is possible for all cases» 88 will be shown in 
following paragraphs. However, one should realize that the portion of 


tangent that osuld be used as an approxi mation to the true adiabatic 


relation is that portion which lies in the first quadrant. Thus the 


upper limit of velocity for practical application of the theory is 


when 3 =0 . By using eqs. (17) and (18), this upper limit 
is found to be 


(9) 


Since the point Ar, 5 ) lies on the true adiabatio ourve, the 
relation hi z b can be used, and eq. (9) becomes: 
D 


SE f 2 4 2 
(Bag SET o 
This relation is plotted in Fig. 2, Since for most practical cases 
it is not likely that the ratio EF) will rise to values much 
higher than 7 , 3 will remain positive, and this theory will give 
an approximate solution, 
If the flow is irrotational, there exists a velocity 


potential f such that 


SE OS 7” 1) 


where u V are the components of W in Y and az 
respectively. The equation of continuity, 


AE) BLE) 


will be satisfied, if the stream function 4 is introduced such 


` that 


N S y= av PR y _ 2f 
d % D 5 2X T 


Now if the angle of inclination of the velocity W to 
the y axis is P » eqs. (11) and (12) give: 


Ab = wr coop di ar m 
dý- W ged dy + wo cug v4 (13) 
Solving for Se and 4 ! 
d: ZA: er ^j " 
y- Muß. d + tip 


So long as the correspondence between the physical plane and hodograph 
plane is one to one, or mathematically 


express y and * as funotions of Aur and so 


as funotions of A . Thus, 


A. di dut ds 


= » di / 
dY = Yor da + dh 
where primes indicate the derivative with respect to variables 


indicated as subscripts. Now substitutingeg. (15) into eq. (14), one has: 


Mi (Ede LEI da HGH, Sg 
dr (Ah 12H Ue IRE TR 


Since the left-hand side of eqs.(16) are exact differentials, one can 


apply the reciprority relation, and obtain: 


PLAE- ER): ee ii 
REH SE 


rn —] - 


Carrying out these differentiations and simplifying with 
the aid of eq. (7), eq. (17) gives: 
== Um — 24 = — e d deu, 5 ; 
Ye ZEN df, m 
$ 
e 


Ge, YP ay le 


Sj 
Solving for pl and 5 A 


d S DES / (19) 


Eq. (19) oan be further simplified by introducing & 
new variable c)  , such that 


fees ae (20) 
Then eq. (19) becomes: A d 
4 7-5 


d o (21) 
4^ fe 
This can be easily recognized as the Riemann-Cauchy's differential 
equations and thus A tar) d must be an analytio funotion of 
-i . However, for convenience of calculation, another new 
set of independent variables A = V/A i d 


VW swf are introduced where W = 4, e” . 


Then eq. (21) can be written ası 


oun xtv) 

ee (22) 
AV ol 

and also by integrating eq. (20), 


W = — a ; > (25) 


ST 


ams 
— E 


a 


ana diri dn W (24) d ' 
4a, W? EN! 
Substituting into eq. (8), we have: Be! 
ha AW Ra | 

$ tal-M (25) Ba 
Eqs.(22), (23), (24) and (25) are the basic equations of 1 A | 
the present theory. One recognizes eqs.(22) as the Riemann-Cauchy A ^ A 
differential equations, and thus the oomplex potential f EP 4 nf E 
must be an analytic function of VA 7 H-V ‚or l 


prip- Fl) HE) 7 H 
= - 26 
bill- En) = En) ? 
Now, it is necessary to find the values of Y and 
corresponds to a given set of values of A and y 9 1.85, 


to find the transformation from hodograph plane to physical plane. 


By using eqs. (24) and (25), eq. (14) oan be written as: 


pa La Ye m 4) 


DU zs 


(27) 


4- MEI, MEE 


E 


2 g. 
where W = H e V . These equations oan be oombined into 


one equation by means of eq. (26). Thus, 


d E Ld 


ha, (28) 


Hence, if an analytic function Hw) is given for each value of 
W , the corresponding real velocity oi oan be calculated by 


eq. (24). Then the coordinate of the point in the physical plane at 


which this velocity occurs can be oaloulated by means of integrating 
eq. (28). The pressure at this point is then given by eq. (2). However, 
there is an objection to this procedure, That is, the investigator 
has no way of telling whether with the chosen fumotion Fi W) » he 
is going to obtain the desired shape of solid boundary and flow pattern 
which he is interested in, In other words, this procedure still suffers 
the difficulty of boundary conditions as it is common to all hodograph 
methods. 

However, due to the particular simple relation of eq. (28), 
one oan ascertain approximately the resulting shape of the body by 


starting with the particular function, 
Hw) = dii y- VA SCH (29) 


where M, is the transformed undisturbed velocity to be interpreted 
by eq. (23), and e 1s the complex coordinate pu . This 
function is so chosen as to give the flow of incompressible fluid over 
the desired body shape in coordinates t and p . Now interpret 
the real velooity in Ss H plane of the incompressible fluid as the 
transformed velocity MW in the hodograph plane for compressible 


fluid, Then it is known that 


W- 4 


= VW de (30) 
Thus = AGE) 
nom Ae (31) 


Substituting eqs. (30) and (31) into eq. (28), one has: 


PE ER ue 


Integrating, 


ME EN 
72477 (32) 


Therefore, the complex coordinate in the physical plane 
of compressible fluid is equal to the corresponding complex coordinate 
in the physical plane of incompressible fluid plus a correction term. 
Since this oorrection term is usually small, the resulting shape of the 
body will be quite similar to the one in incompressible fluid. The 
factor in the correction term depends upon the Mach's number of the 
undisturbed flow only. This oan be seen by means of eqs. (7), (8), 


and (23), because from those equations one has the following relation: 


EC NEN. NE 
Zu Dm (33) 


where a is the Mach's number of the undisturbed flow, 
To calculate the velocity in physical plane, W is 


first obtained from eq. (30) and then with eq. (23): 


Fa) UI — mi 
Frl (34) 
77m, pp IW JE 
FEB 
Then the pressure can be caloulated by using eq. (2). 


With some calculation the following expression for pressure at any 


point is obtained: 


E Erf res" (36) 


ag 


In this part, the theory developed in Part I will be 


applied to calculate the flow over an elliptic cylinder at zero angle 


of attack. The incompressible flow over an elliptic cylinder in the 
complex coordinate can be obtained by applying the Joukowsky's 
transformation to the flow over a circular cylinder in the complex 


coordinate z' with its center located at the origin. Therefore, 


the function FIN) or WG) can be written as: 


(37) 


m! It is convenient to carry out the computation by using the Z i 


coordinates. Thus eq. (32) is rewritten in the following form: 


H ge Dë EE WHEE T 
Aa 


LN If one limits one's self to caloulate the conditions over NW ! l 


E the surface of the elliptic cylinder only, 
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where A is the argument as shom in Fig. 3 and Z is the radius of 
the circular section in the g= plane which determines the thickness 
ratio of the elliptic section in the Z-plane, Substituting eqs. (36)» 
(37) and (39) into eq. (38), and carrying out the integration, the 
following expressions for the X me f coordinates corresponding to 

e are obtained by separating the real and imaginary parts, 


= M: "VE, 
LEO ob - 5 Me IM boat ru eA 


(40) 


Jn [i Zei Je. f bt) wid HED * HS. 


The horisontal and vertical semi-axis of the near elliptic section oan 
then be oaloulated by substituting £ =p and f= Fs 


respeotively into eq. (40). The thickness ratio f is thus obtained 


2 p= (fig) eie ae) 4a, A Co, AC 
"Ire MONS ated] 


For a given thickness ratio and Mach's number of undisturbed stream, 


2 
the value of $F is first computed by means of eq. (33), and 
then eq. (41) is solved graphically for A 


(41) 


After f is obtained, the coordinate A and Pj for 
each value of / can be computed by using eq. (40). It is fortunate 
that the values of A, so obtained lies very close to the true 
elliptic section as shown in Fig, 4 Hence, the velocity and the 


pressure distribution obtained by using eqs. (34) and (35) are considered 


x j - | es erm. 4 se F: RE 4 Vf 


" 


0437 


dics sal 


e 


Ay 
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as those over the true elliptic sections. 

Calculations for two thiokness ratios, C= 05 
and {= 0. » are carried out and the results shown in Fig. 6, 
Fig. 6, Fig. 7 and Fig. 8, together with those of Kaplan (Ref. 4) and 
Hooker (Ref, 10) inoluded, In order to oompare the results with that 
obtained from less exact theory of Glauert and Prandtl, these compu- 
tations are made (Ref, l ), and the results also inoluded in Fig. 5 
to Fig, 8, It is sem from those figures that the effeot of compressibility 
of fluids shown by the present theory is more pronounced than those 
shown by Kaplan's and Hooker's caloulatim. In view of the fact that 
Kaplan and Hooker carried their oaloulationmly to first approximation, 
and higher approximations will inorease the value somewhat, it is 
probable that the values given by the present theory are nearer to the 
true value. As to the results given by Prandtl-Glauert theory, it is 
seen that the theory breaks down near the stagnation point and the 
results are poorer for thick sections than those for thin sections, 
However, in oase of thin sections this very approximate theory does 
give satisfactory information for a point not too close to the stagnation 
point, 

If the velocity of flow over a body is gradually inoreased, 
the maximum local velocity in the field will also be increased, hen 
a certain velocity of the undisturbed stream is reached the maximum 
local velocity reaches the local velocity of sound. Then shock waves 
usually appear and the drag of the body suddenly increases. Therefore, 
this velocity is of considerable interest to practical engineers and 
is usually referred to as the critical velooity of the body. It is 
shown by Kaplan ( Ref. 4) and others that at this oritical condition 


the ratio of maximum velocity Wy in the field to the undisturbod 


velooity W, , is related to the Mach's number 
stream in the following manner: 


2 


Kay in the flow over an elliptic oylinder at zero angle of 


attack ooours at the top of the cylinder. Using eqs. (34) and (33) 


r A4, 
LA the value of en is found to be: 
/ 


l a 
7 pi f Li- 827] a" (43) 
AB n HF UEM än 


Equating eqs. (42) and (43), the equation for caloulating 


40, 
the critioal Mach's number E) of the undisturbed stream for 
sach value of L is: 
sh” 2 
STH 


£4.44 — nn — 
PELHA LIU 


A n f AL 27773 (44) 


" 
f th 
Xnowing Cor for each value o A ‚the 
corresponding value of d can be oaloulated by means of eqs. (34) 
and (41). Fig. 9 shows the result of this oaloulation with Kaplan's 


da ee 


or e (49) 


4 


Therefore, the quantities R and A have the dimension of a 
velooity and so R and > can be considered as a new velocity 
in the plane whose coordinates are denoted by YY ana V . This 
relation between the y 2 plane and the P4 taċ y plane is shown 
in Fig. lO. It is thus sem that if the undisturbed stream in the 


D A4 plane is in the positive y-direction, the undisturbed stream 


in the 4 * plane will be in the negative X direction. Furthermore, 


if Ü 
20 ER = r 
5. p us | 
M pt rer? (50) 


one has from eq. (48), 


(51) 


* x the LY plane, as shown by H. Bateman (Ref. 11). 


(db 
r ago RE PN ies o etc EEE m 
Comparing qe (51) with | eq. (12), it is evident that 2 
H 
^ T? 
ei can be considered as the density of a fluid in tho new ERS plane, MA 
| l So far the relations obtained are general, i.c., they 


| &pply to fluids of arbitrary properties. However, sinoe only the flow 
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z of incompressible fluid has been studied exhaustively, it would be CH 
+ interesting to find out what will be the flow in the ZE plane if a 
| ! 5 the flow in the Y 4/ plane is incompressible, If the fluid in the e 
he kl Wi plane is incompressible, then f-/ » and the Bernoulli's 
head is: 


wa Bau’) 
ji (52) 


2 
Now in the VY plane, let f denote the pressure and Q denote , 
£^ $ 2 3 then Bernoulli's theorem gives: | 
' 


| : 
SE r4 vgs constant (53) 


In view of eqs. (49), (50), and (52), eq. (53) can be 


written in the following form: 


/ MG) AP hl 
EEE 


By differentiating eq. (54) respeot to iG 3 


4 = constant (54) 


N multiplying the Ke expression by [3 and then integrating 


with respeot to = » one has the Robe relation connecting 


the pressure P and density 0 for the fluid in the PA V plene: 


P = MALA SCH Zo = (55) 


Comparing eq. (55) with the approximate adiabatic relation 


eq. (2), also noting eq. (16), it is evident that eqs. (55) and (2) 


are identical if 
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th at p Ze 
it Aor Ul (56) 
and Constant = fa Zë VA a 


In eq. (56) 4 is the velocity of sound in tne P H plane, and the 
subscript 1 refers to the conditions in the undisturbed parallel 
stream, Hence, 4. is the Mach's number of the undisturbed 
parallel stream. i 

To find the coordinates E 4 and Lé in terms of 
X, » eq. (48) has to be integrated, It is convenient in this case 
to use the complex potential of incompressible flow in the L 2 plane, 
Thus, if 


prif- w Glatg)= ar GU) 


(57) 


Then it can be easily shown with the aid of eq. (52) that 


deli! dae 227. (58) 


It is interesting to notioe that eq. (58) is identical to the 
correction terms in eq. (32) except a factor. However, there is a 
fundamental difference between the method explained in Part I and the 
method explained in this Part. The previous method will give different 
solid boundaries for different Maoh's number of the undisturbed flow 
even if the same function (7 is used. Eq. (58) shows that if the 


kd Zn A , H " , PL GE 2 


pnm el 
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A same function G is used, the resulting solid boundaries are all i be 
n geometrically similar. 1 n 
V au By using eqs. (52) and (49), the components of velocity Y: | 
eU in the LY plane can be EEN m E 
Lo f- E bi 
1 


E " TRE B 


c P lija tee) 


i 
$0, 1 

= / H LE n a 
wi Bei E Af 

$-e Lip e» 

ak i aC N 


In eqs. (59) and (60), it is assumed that the undisturbed flow in 
the kä plane is in the Y —direction and so tne undisturbed flow 
in the pa y plane is in the K direction. The relation between 


Ai: and A can then be obtained from eqs, (56) and (57), that is: 


tien EM Pa - 48) 


ER 


Ke 


(61) 


Thus, eq. (60) oan be rewritten as: 


- Wr T " BL 7—- (62) 


For a given condition of flow in the y Y plene, eqs. (56) and (61) 


will determine the characteristios of the flow in the L plene. 


Ca Using eqs. (55), (56) and (50), the pressure in the pA y plane oan be 


evaluated as: 


` = JF hel 

T » L 

| : (63) DA 

ki $ , 
|n i i 
| Thus as the method is developed one starts with a certain funotion Wb 


G2) representing the incompressible flow in the KE plane, we 
| first calculate the coordinates ZZ r and thus the solid boundary 

l in the P4 d plane. Then with the desired value of b) for flow 
in the Dé Y plane, the velocity and pressure een can be 
calculated by means of eqs. (61), (62) and (63). Thus the flow problem 


of approximately adiabatic compressible fluid is solved by means of 


an auxiliary function 6 (2) for flow of incompressible fluids. 


| 
However, this procedure suffers the defect of difficulty in obtaining 
| the desired solid boundary compared with the first method developed A 
| in Part I, because now SC? is not related to Z bya small correction m 
as is in oq. (32), but by a transformation given by eq. (68). t 
3 However, a much more useful form of eqs. (62) and (63) 
l can be obtaincd if one remembers that the solid boundarv of the flow t 
in the Y v plane is independent of Mach's number, but only on 


Gle) » or in other words, on « Now substitute eq. (61) 
A 


into eq. (62) and let LE) —o 
Ba Ki 


or 


od 
where m, denotes ‚the velocity for flow in the 24 r plane 


when (£) =D . Then substitute eq. (64) into eq. (62). The 
following relation between the velocity of incompressible flow and 
the velooity of compressible flow over the same solid boundary is 
obtained: 


£)- t£] d 


(65) 


where E? is a function of fi) given by eq. (61). Similarly 
eq. (63) can be written as: 


Q E 
= (177) at (66) 
Gd 


f Jf for fi 4 — is denoted by T, , then it is well 
7 


known that 
l-I 
or J=,- 7, (67) 


Substituting eq. (67) into eq. (66), we have: 
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Eqs. (65) and (68) together with eq. (61) thus enable 


one to calculate the velocity and pressure distribution over a body 
in compressible fluid once the velocity and pressure distribution over 
the same body in inoompressible fluid is mown. The latter can be 
obtained by either well-known analytical means or by experiments. 

If eq. (61) is substituted into eq. (68) and only first order 
terms of LE) 1s retained, then 


UA $ 
adi B $ 


Thus, if e and CG are the lift coefficients of an airfoil T 


o (€) and at HI =0 and the same angle of Fu 


attack respectively, one has from eq. (69): 


o 
} I mi 
3 W 1-4) 


a = which checks with the Prandtl-Glauert theory (Ref. 1). 


f B | 
| 

By inverting eq. (68), one has: | de 
k 

"i 


N — EM 
r pA ZE aw (71) 


If the theory developed in Part III is true, then the value of VA 


calculated from Z4 at different values of Mach's number fi Si of 
7 


k ' undisturbed flow should fall on the same ourve, This is done for the 


qe ee 
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case of upper surface of NACA 4412 airfoil tested by J. Stack (Ref. 2) 
at -09/5 á angle of attack. The 7] are replotted in Fig. 10, while 
the FA curves calculated from y/4 curves are shom in Fig. 11. 
The derived values of Kë Lie fairly olose to me ourve. Since the 
data for this calculation is taken from a small scale figure, the 
deviations must be largely due to inaccuracies involved. Fig. 11 can 
thus be taken as an experimental check of the theory. 

One of the immediate possible applications of the theory 


developed in Part III is to determine the critical speed of a body. 


H l Sch and Z, min “re the maximum speed ee 


C minimum pressure over a body. Then by eqs. (65), 


A ume 


(87) and (42), the equation for eritioal Mach's number is obtained as: 


(e - TEL FH o» 
l | Ft Fol 72 
n laly ^ tl, HT 


Wei 


or 


Te 3 


6, 
In both eqs. (72) and (73), f is a function of Fl given 


by eq. (61). Eq. (73) is plotted in Fig. 12, together with a curve 

given by E. N. Jacobs (Ref. 13) based Ny EAN, theory. It 
4s thus seen that Jacobs's ourve sives a oritical Mach's number, which 
is a little too high. This difference is, of course, to be expected, 


due to the higher order approximation of the present theory. 


2. 


3. 


4. 


6. 
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Approximation to the adiabatio relation by means 
of a tangent. 


Relation between the maximm velocity (at which the pressure 
is zero) and Mach's number, 


Variation of the parameter of transformation from incompressible 
flow to compressible flow with Mach's number, 


Notations used in calculating the flow over an elliptical 
oylinder. 


Flow over an elliptical cylinder with thickness ratio = 0.5 
at Mach's number = 0,5, 


(a) Velocity distribution 


(b) Pressure distribution 
Fig. 6. Flow over an elliptical cylinder with thickness ratio = 0.1 ] 
at Mach's number = 0, 857. 
kI (a) Velocity distribution 
I" (b) Pressure distribution ' 


gi Fig. 7. Variation of critical Mach's number of an elliptical cylinder 
Bn with thickness ratio. 


F Fig. 8. Relation of the velocity components in the plane and 
plane, 
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SIMILARITY LAWS OF HYPHRSONIC FLOWS 
by Hsue-shen Tsien 


| Introduction 

Hypersonic flows are flow fields where the fluid velocity is 
much larger then the velocity of propagation of small disturbances, 
the velocity of sound. Th. von Kárman (Ref, 1) has pointed out that 
in many ways the dynamics of hypersonic flows is similar to Newton's 
corpuscular theory of &erodynsmios, The pressure acting on an in- 
clined surface is thus greater than the free stream pressure by a 
quantity which is approximately proportional to the square of the 
angle of inclination instead of the usual linear law for conventional 
supersonic flows. E, Sanger (Fef. 2) has, in fact, used this concept 
to design the optimum wing and body shapes for hypersonic flight at 
extreme speeds. 

Recently, von Kerman (Ref. 3) has obtained the similarity laws 
for transonic flows where the fluid velocity is very near to the 
velocity of sound. He deduced these laws by using an affine trans- 
formation of the fluid field so that the differential equations of 
the flows are reduced to a single non-dimensional equation. In 
this paper, the same method is used to derive the similarity laws 
for hypersonic flows. These laws will be, perhaps, useful in oorre- 
lating the experimental data to be obtained in the near future by 


hypersonic wind tunnels now under construction. 


^ Differential Equation for Hypersonic Flows 


If u, v are the components of velocity in the x, y directions 


and a is the local velocity of sound, the differential equations for 
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irrotational two-dimensional motion are 
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Now if a slender body is present in an otherwise uniform stream of 


velocity V in the x-direotion, equations (2) is satisfied by intro- 
duoing the disturbance velooity potential f defined as 

a= Vt * 

v= 2% 
If a, is the velocity of sound for the gas at rest and a9 is the 


(3) 


velocity of sound corresponding to the free stream velocity V, 
then there are the following relationst 
ài. .* qr 2 4 £ 2 OY j2 
a (wv) DV ve PEP EI] qa 
ei Ca sud 
a= a- HV 
where A is the ratio of the specific heats. 
For hypersonic flows over a slender body, both a° and 28 VA 
aX, op 
are small in comparison with V. By substituting equations (3) and 
(4) into equation (1) and retaining terms up to second order, one has 
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(8) 
Pr "D A no 


Here M° is the Mach number of the free stream or 


Mz V (6) 


gës? 


Similarity Laws in Two-Dimensional Flow 
If 2b is the length or chord of the body and d the thickness 


of the body, the non-dimensional coordinates 4 and 7 oan be de- 


fined as il 
X-b£ MA. 
A (7) n 
= b > 
where n is the exponent yet to be determined. von Kármán (ef. 1) TT 
un 
has shown that for hypersonio flow over a slender body the variation > 


of fluid velocity due to the presence of the body is limited within 
a narrow region close to the body, the hypersonic boundary layer, 
Therefore, in order to investigate this velocity variation, one must 
expand the coordinate normal to the surface of the body. This is - 
similar to the case of ordinary viscous boundary layer, where Prandtl's 
simplified boundary layer equation is obtained from the exact Javier- 
Stokes equations by a coordinate expension normal to the surface 
of the body. From this reasoning then, n must be positive so that 1 
is much greater than G)/(). This surmise is substantiated by the 
later calculations to be shown presently. 

The appropriate non-dimensional form for the velocity potential 
Y is 

Ọ= ab x A, 7) (8) 


By substituting equations (7) and (8) into equation (6), one has 
17 Me ag LOT alle A aa E i 
ed Seier? - 
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The boundary oonditions at infinity require that the flow 


velocity be V, Thus 


Z- # =) at ^^ (10) 


If the slender body is a diu one, then the oondition at the 


surface of the body oan be written as 


PT ^ 


where n(f ) for -1< < 1 is a given function describing the thick- 


ness distribution along the length of the body. Equation (11) oan 


be converted into the following form by means of equations (7) anà 


(8). 2 tn 
2] m we) d) (12) 


Since the body is thin, Käl is very small. Therefore the first 
group of terms in equation (9) is negligible in comparison with the 
rest, Then both the differential equation and the boundary conditions 
can be made to contain only a single parameter if one sets 


ns L (13) 


ME=K 
then equation (9) becomes 


je Iit - ary BOTH - l a a (14) 


That is, if 


and the boundary conditions become 


pe 3 = at d (15) 
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and 
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The meaning of this similarity law is the following: If a series 
of bodies having the same thicknsss distribution but different thick- 
ness ratios (s/s) are put into flows of different Mach numbers 
MO such that tho products of MO and (5/0 remain constant and equal 
to K, then the flow patterns are similar in the sense that they are 
governed by the seme function //§, 7) » determined by equations (14) 
and (15), 

If po is the stagnation pressure, p9 the free stream pressure, 


and p the local pressure, then 


ĦA a n] 
sey du 
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In notations introduced previously and by retaining terms of proper 


m&gnitudes, one oan write the expression for tne loosl pressure as 


PT uo B £2] 77 
$ A [1 ie RS) y? (16) 
The drag D of the body can then be caloulated. It is given by the 


following a 
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If one wishes to compute the drag coefficient Cp, then one 
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For a given thickness distribution, the quantity within the brackets 
is only a function of K, the similarity parameter. Therefore, one 


can write 


c,- 55 A(K) ~ Alm) 


Similarly, one obtains for the lift coefficient D of the lift L 


(18) 


the following law: 


L , Ye SAME (19) 
C = EV) D A ) HAM 5 ] 


These similarity laws show that for bodies of the same thickness 


distribution at angles of attack proportional to the thickness ratio 
( Wé ), the quantities (cgi?) and cue?) are functions of the 
single parameter K * M9 £ A 

Equations (18) and (19) agree with the results of the more 
limited linearized theory of Ackeret (Hef. 4), Aooording to this 
theory, for similar bodies in the sense stated above the drag 
coefficient and the lift coefficient are given by 

CECR 


D 


Lei 
C, A o E. 
Mr 
For hypersonio flows of very large values of M?, these expressions 


e HK 


reduce to 


(20) 


b (21) 
C, ~ Lg 5 


Equations (20) and (21) agree with equations (18) and (19). Equa- 
tions (18) and (19) are, however, more general and complete. 


Axial trical Flows 

For axially symmetrical flows, the ordinate y is the radial 
‘distance from the axis to the point concerned, Then a similar 
analysis leads to the following differential equation and boundary 


Ji conditions! 
Dir \ 
l M mi Ly i : 244 
ia a e AH 
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] ] H ke 
l 
Ke (24) 
d DA), 3 Kt) É Se 
d 1 j where a(g ) is the distribution funetion for oross-seotional areas 
im t along the length of the body and K » MO i as for two-dimensional 
N flow 
l The drag coefficient Cp referred to the maximum oross section i 
ki of the body is then governed by the following similarity law: 
Y C = us AUT) (26) EI 
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The "Limiting Line" in Mixed 


Subsonic and Supersonic Flows 


of Compressible Fluids 
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It is well-known that the vorticity for any fluid element is constant 


if the fluid is non-viscous and the change of statew of the fluid is isentropic. 
When a solid body is placed in a uniform stream, the flow b pen of the 
body is irrotational, Then if the flow is further assumed to be isentropic, 
the vorticity will WR eee zero over the whole field of flow. In other 
words, the flow is irrotational. For such flow over a solid body ef-firite 
dimesretzrm, it is shown by T. Theodorsen (Ref. 1) that the solid body exper- 
iences no resistance. If the fluid has a small viscosity, its effect will 
be limited in the boundary layer over the solid body &nd the body will have 
a drag due to the skin friction. This type of essentially isentropic 
irrotational flow is generally observed for a stream-lined body placed in 
a uniform stream, if the velocity of the streamiis kept below the so-called 
"oritical ef? 
aż At the "critical speed" or rather at a certain value of the ratio of ZA "oo 
Aft : a" 
We Welocity of sound, shock wavés appear. This phenomenwm is called the 
"compressibility bubble", Along a shock wave, the change of states of the 
fluid is no more isentropic, although still adiabatic. Thés results in an 
increase in entropy of the fluid and generally introduces vorticity in an 
ċrżġinaliv irrotational flow. The increase in entropy of the fluid is, of 
course, the omsequence of changing part of the mechanical energy into heat 
energy. In other worüs, the part of fluid frrected by the shock wave has 
much less mechanical energy. Therefore, with the appearance of shock waves, 
the wake of the stream line body is very much getigert, and the drag increases 
drastically. Furthermore, the accompanied ehange in the pressure distribution 
over the body RR Dd rm abba ta it and in case of an airfoil 


decreases the lift force, 


of extremely large viscofs effect. In other words, the boundary layer 
will be so thick as to occupy the main portion ef the nosvle pessage. To 


demonstrate tnis effect, let us-sensider—that the length of the test section) A 
leta L and the width of the square test section be b. Then the Reynolds num- 


ber based uon the conditions in the test sestion is Re = D where U is 
the velecity in the test section. If, as a rough estimate, we take the 
thickness of the boundary layer to be sero at the beginning of the test 


seotion and equal to a value $ caleulated by the well-known Blasius formula 


for & flat plate at the end of the test section, then 


= L- 
How if this boundary layer «ctually oocu;ien half the tumnel width b/2, 
then 
>. L 
KA As (2) 


Qu the other hand, tbe ratio of the mean free path Í amd the boundary 
layer thickness } is known (Ref. 2) to be equal to 


L- Kr ed Wi 


where 7 is the retio of speeific heats and cam be taken as 1.4, M is the 
Mach mmber in the test section. By combining (2) and (3), va have 
ote ) Vie 0.0557 M (4) 
c.g tt /This relation is shown in Fig. 1. Thus for a Mach mmber M equal to 2, 
äi f and L/b = 2, the boundsry layer will completely fill up the test section, 
7. 4f the mean free path is equal to 3.68 of the Boundary layer thickness or 
2.9% of the tunnel width. This means that the extremely strong viscous 
effect at lew derivities makes the ordinary soncest of designing 3 wind 
tunnel totelly inapplicable, 
Tbe tene Gest con be also demonstrated by calculeting the ratio of 
friStiomal loss on the walls ef the test section amd the shock loss in the 
diffuser after the test sectien. Consider the diffuser to be a straight 
tube of approximately the same crese-sestional ares as the teat section, 
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Introducing the mean free path ratio given by (3), we have 


A 6.52894 


t = 


d: z (8) 
> CO Te] 
This relation is plotted in Fig. 2. Therefore if the Mach number M is 2, 
and L/b = 2 as before, then when the ratio(t/j) is 0.056, the ratio of 
frictional loss to shock loss is 0.628. Thesefer® the frictional loss and 
>the shock loss is of the same order of magnitude, 
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These elementi:v colculationsYmekes it clear thet for the design of the 


nozzle and test section for à superserodynamics wind tunnel, it is ao longer 
possible to separste the compressibility effects and the viscous effect. In 
fact, the concept of boundary layer is also of doubtful value due to the 
extremely small Reynolds number encountered. Therefore to actuelly design 
such a nozzle to obtain the nearest approximation to the idesl uniform flow 
it will be necessary to use the exact MNevier-Stokes e.ustions insted of 
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wdd:tlonal corrections murt be :dded (Ref. 2). however, recent investi a- 


d 

"n tne up roxizzte bound:.rv layer ecu: tions. Of course, it muy be «argued that 
P for Guer erodynanics, the tavier-Atokes e uati os for no more exact end 

l 


tions by R. ċchumberg (Ref. 3) have shown that these addition=1 corrections 
are smell in case of slip-flows concerned here and will not essentially 
L alter the flow pattern, Hence for a first approximation just like the non- 
l viscous isentro;ic flow as a first approximation for ordinary supersonic 
l uoisles, we can use the Navier-Stokes equations. The simplest case to be 
considered is certainly the axially symmetrie noszles. If x is the coordinte 
in the axial direction, Y the coordinste in the redial direction end u 


and v sre the corresponding velocity components, ( H the equations are: 
Hy LE) =o | (9) 

gt = - 2b EH (20) 

poy -= M + Grad (T), (2) 


SH + er) -8- JuGradlc), t 0 Grad lr) d AE. trir zv (12) 


where T eS s BA 
plo a (Er 
& = dessipition function 
© T = stresses tensor 


(9), (10), (11) and (12) togetner with the equstion of states 
4 =R? (23) 


then determines the five un«nows u, Y, faf , and f. Of course, the actual 
process of making this celoulation will be extremely tedious and some a; rox- 
imation sethod oi solution may have to be developed. One possibility would 

be to adopt the Karman-Pglhausen method for boundary lsyer to this ccse: We 
intkgrate the differential e wations ence with respect to r and thus only 

try to satisfy the equations "on the Xveraget over the cross-section of the 
nozsle. The "distribution" of u, v over the cross-section will then be set 
in the form of c colvnomisal inr. Initial study in this procedure is alresdy 
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Bude by 5. A. Schaff (Ref. 4) at the suggestion of the author. 

For orcinury supersonic diffuser, high efficiency of pressure recovery 
can be generally .chieved by using a long diffuser. However, for surersero- 
dyn.mice wind tunnel, due to the extremely l=rge loss through friction, long 
diffusers sre undesirable. In fact, the pressure lose c-n be reduced by 
using r shortest ;ossible diffuser. 


jl 

2. Flow Messurement | 
The .uantities which determine the flow field ars t«o out of the | 

three varisblesf »® , T and the velocity com,onente. The ,uantivies , Pe y 


T sre related by tne equstion of states snd tierefore only two is necessary 

for the detersination of all three, Generaliy for wind tunnel work, the 

quantities sctusliy meueured ere f, f and b the magnitude of the velocity. 1 
zor the me:surement of preseure, a afsnomet r is used. For ordinary 

pressure, one uses a fluid z nemeter filled with water, alcohol or mercury. 

However for tne extremely low pressure encountered in the WsperKérodyc mE 


flow, some other form of manoweter is necessary. One of the most weed'type 


| is the Pirani gatge. The conventional form of Pireni ygsages hes + pressure 
| sensitivity of about 10 micron,® It utilizes the eh nge of temper«ture 
| of u wire hected with constant energy coused by a clu.nge in the pressure of 


the gus surrounding it. The temperature chenge 1s measured by the ch nge in ġ 
| 
| 


the resistance of the wire. The wire is locuted in 5 smali ch mber which 
is connected to the „oint oi measurement by a hole, flush with the gas stream 
if static pressure is to be measured. The question of best design of the 
| connecting tube for uuick response is studied by S. 4. Schaft. (Ref. 5). 


To mecsure the density ¢ » the conventions! method utilises the differ- 
ence in the velocity of light rays in mediums of different 49DBÍty- ich 
| different optical errengements, we have the shadowgraph method, the schlieren 


method :nd tie Interferometer method. However, if the density of the "9dlum 
1s very low ss the case of superaerodynamie flow, the sensitivity of these 


methods become extremely poor. For instence in c.se of schlieren method, the 
percentage change in illuminztion I by psssing through : region of thicknese 


D 
* 000 micron = l m m. Hg. one abes. = 2b € i een 


b 1s given by 


AI > 
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where WI the air density at 32°F and 1 «tases. here ;ressure, and 4P/an 

is the density gradient normal to the light ray. f and zre the local 
length and the normsl unobscured width of the light source image . erpendic- 
ular to the suife edge. k is factor of order 1, determined by the particular 
o: tical ‚atn used. herefore the sensitivity of the schlieren method de- 
creases with the factor ( WA ). Some imporvement c:n be made by altering 
the quantities f sud E but practical limitations and diffraction difficulties 
do not tilow tie incresse of sensitivity to sutisfactory values. 


(14) 


A new upproach to this ;roblem density meunsurement is the method of 
»bsorption. It is found for instenoe, that oxygen at low pressures shows 
a strong absorption band at wave lengths around Kaz n or ultra-violet 
light. The percentage :bsorption is proportional to tie number of molecules 
that meet the light ray and is, therefore, proportional to the density of 
the gas. The measurement is then similar to that of tne interferometer 
method where the density of gas is determined. A similar method is the 
utilis-tion of tne after-flow of nitrogen. These methous «re now being 
studied by R. A. Evans (Ref. 6). 

The conventional method for the meusurement of velocity is through the 
use of dynamic pressure rise in a ritot-tube. A straight forw:rd ap; lication 
of this method is, however, difficult for reretíed gases. The formule used 
is, however, b.sed upon the neglection of viscosity effects. But for rare- 
fied gases the viscosity effect im of greet importence rg ‚ointed out in 
the previous section. Then the dynamic pressure would be ,uite different 
than that given by the usual formula. To estimate this effect, let us 
consider the case of low Mach number so thet compressibility effects can be 
neglected. Then us à first approximation, take the flow field eround the 
Pitot-tube as tist of a source of strength $ in non-viscous flow of uniform 
velocity U. (BEE). The "radius" of the tube a is 


end the stegnetion ‚oint is located at 


Mur 


R NC = 4 , (15) 


The velocity introduced by tne source is then 


By celcul:tiag the viscous stress from this spproximate disturbence velocity, 
we tive for flow slong the axis 
JE xl Via tS tL (16) 


t? De 
Hence if b is the u ‚ressure end A the static pressure, 
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EE 


Or Sm E I (17) 
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For rsrefieċ gŁses, the value of Y/al/ or the reciprocel of the Reynold's 


nuaber of the Pitot-tube could be of the oraer of unity. Then the dynamic 
pressure rise A -F is not the usual values dev but & velue much less 


than that. V Ve-sheli— be sepieusiy in-error-i-ve-use-the-erdinary-fornula 
he-ealculate the-welooity-7/ . 

when the velocity of flow is high, we have the added com licetion due 
to tne shock. Tne conventional Rayleigh formua for Pitot tubes in super- 


sonic Ilow is based upon the usżumption of very thin shock weve checd of 
the Pitot tube. Now the thickness of the shock ia proportionti to the mean 
free path. Hence in rarefied flows, the thickness of the shock will be so 
increaset &5 to cause interference with flow in the neighborhood of the 
Pitot-tube. This together with the viscous effect zentioned in previous 
paragraph definitely show the inapplicability of the Rayleigh 1ormul: for 
DEE velocity of rarefied gases, 
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bet measuring device to supercerodynamic flows, one is naturelly let % the 


thought of owner avenues of ap;roach, One possibilit; is the use of ihote 
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j Parameters of Flow 

The two parameters thet are directly connected with the flow field 
are the Reynold's number Re, defined as 

ie = d 

where y'is the xiuemrtic viscosity, L the typicel lineer disension of the 
body; end tne Mach number Dy of the free strean. This is true even for 
slip flows and free molecule flows due to the fact thet the ratio of mesn 
free path to the typical dimension can be also expressed in terms of the 
Reynold's number «nd the Mach number, 

However, as the pressure or density is reduced, the solid boundery of 
the flow enters ectively into tne flow conditions by re.uiring not only that 
the microscopic stream velocity be t*ngential to the surface but that the 
interaction of the molecules end the wall be considered end thet the radis- 
tion of energy to and from the wall be taken into account. The interaction 
of the molecules with the wall is so far exnressed through the fraction 5 
of molecules thet are diffusely re-emitted from th rand the secoumedation 
coefficient X . (We have reasons to believe Laf {and Y are functions 
of the tez.erature of the wall aad KE ty distribution. 
Therefore the interaction of tne molecules with the wall is the same only 
if the wall teupereture, the ges tem;ersture and the Mach number of the 
gzs above the wall is the same. These considerations seem to indicete then 
that ‘or the model test to be similar to the prototype, the model rust be 
made of pame surface material as the prototyje, and the fluid must be the 
Bane, onu jurtiermore the foliowing parameters must be the same: 

(i) Reynolds number Re 
(2) Mech number M° 
(3) free stream temperature, T° 
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The r:di.tlon hest loss from the surface is equal to /6/, per unit 
area. However, if the model is surrounded by the wella of the test chember, 
there is sls0 an bet input due to redi. tion from walls of the test section 
to the model. Let us call this quantity 2. , Then the net heit loss cer 
unit «res of the surfsce of the model is £67 - f , This quantity can be 


D 


rendered non-dimensional by dividing it by pu (bp 7°) >» Call t is new 


porameter /„ , then 
oe et, kj 34 
PU eT) Im (85) 


For the prototy:;e, the heat from the walls of the test chamber is ubsent 
but there may be solar rediation and the radiation from the earth and sur- 
rounding etmos here. (Ref. 9). Denoting tnis amount by d » then the 
parameter A for the prototype is 
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In order for the flow to be same also with res.ect to the raülation heat 
transfer, 


A^ 7A. (23) 


Because of the previous conditions on the Reynolds number und free stream 
temperature, (27) is the same as 
l Burn da p 
fs er he e PM 
where M is the trpicsl Linger dimension of the model and L is the typical 
linear dimension of the prototype. This mecns thet the wall temperature of 
the test chamber aust be so controlled that f, s:tisfies (28). 
This set of rather strict simileriv rules for model testing in surer- 
serodynemic flows is certainly difficult to satisfy. In what way the rules 
can be relaxed is the problem of future rese:reħ. 
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THE "LIMITING LINE" IN MIXED SUBSONIC AND SUPERSONIC FLOW 
“OF COMPRESSIBLE FLUIDS SSS 


Hsueeshen Tsien 
California Institute of Teehnology 

It is wellekmgwn that the vorticity for any fluid element is constant 
if the fluid is non-viscous and the change of state of the fluid is isentropic. 
When a solid body is placed in a uniform stream, the flow far ahead of the 
body is irrotational. Then if the flow is further assumed to be isentropic, 
the vorticity will be z ro over the whole field of flow. In other words, 
the flow is irrotational, For such flow over a solid body, it is shown 
by T. Theodorsen (Ref, 1) that the solid body experiences no resistance. If 
the fluid has a small viscosity, its effect will be limited in the boundary 
layer over the solid body and the body will have a drag due to the skin frio- 
tion. This type of essentially isentropic irrotational flow is generally 
observed for a streem-lined body placed in a uniform strem, if the velooity 
of the stream is kept below the so-called "oritical speed", 

At the "critical speed" or rather at a certain value of the ratio of 
the velooity the undisturbed flow and the corresponding velooity of sound, 
shock waves appear, This phenomenon is called the "compressibility burble". 
Along a shock wave, the change of state of the fluid is no more isentropio, 
although still adiabatic. This results in an inorease in entropy of the fluid 
and generally introduces vorticity in an originally irrotatioral flow. The 
increase in entropy of the fluid.is, of course, the consequence of changing 
part of the mechanical energy into heat energy. In other words, the part of 
fluid affected by the shook wave has much less mechanical energy. Therefore, 
with the appearance of shock waves, the wake of the stream line body is very 


much widened, and the drag increases drastically. Furthermore, the accompanied 


change in the pressure distribution over the body changes the acrodynande T 


moment acting on it and in oase of an airfoil decreases the lift force. T. "n 

All these consequences of the breakdown of isentropic irrotational flow DM 
are ponerally undesirable in on Lied aerodvramios, Its occurrence should be 
delayed as much as possible by modifying the shape or contour of the body. 
However, such endeavor wil! be very much facilitated if the cause r the 
criterion for the breakdown oan be found first, 

Criterion for the Breakdown 
of leentropio Irrotational Flow 

G, I, Taylor and C. F, Sharman (Ref, 2) calculated the successive 
approximations to the flow around an airfoil by means of an electrolyte tank, 
They found that wha the maximum velocity in the flow reaches the local vel- 
ooity of sound, the convergence of the successive steps seems to break dom. 
This fact led to the identification of oritioal speed or oritical Mach number 
with the Mach number of the undisturbed flow for which the local velocity at 
some point reaches the local velocity of sound. However, there is no mthe- 
matical proof for the coincidence of the critical Mach number so defined and 
the breakdown of isentropic irrotational flow. Furthermore, such a definition 
for oritioal “ach number implies that a transition from a velocity less than 
that of sound, or subsonio velocity, to a velocity greater than that of sound, 
or supersonic velocity, does not occur in isentropic irrotational flow, On 
the other hand, Taylor (Ref. 3) and others found solutions for which such a 
transition occurs. Furthermore, A. M, Binnie and S. G. Hooker Lef, 4) have 
shown that at least fo ips case of spiral flow the method of succe sive 


approximation is a convergent one even for supersonic velocities. With these 


facts in mind, one can conclude that the identification of critical speed 


| 
| 


m 


with local supersonic velocity cannot be correct. 

Taylor's investigation on the spiral flow (Ref. 3) indicates that there 
is a line in the flow field where the maximm velocity is reached and beyond 
which the flow cannot continue, W. Tollmien in a subsequent paper (Ref. 5) 
called such lines as the “limiting lines". The velocity at the limiting 
line is never subsonic. However, the true characteristics of such limiting 
lines and their significance were not investigated by Tollmien at that time. 
Recently, F., Ringleb (Ref. 6) obtained another partioular solution of isen- 
tropic irrotational flow in which the maximum velocity reached is approximately 
twice the local sound velocity. For this flow also, a limiting lino appeared 
beyond which the flow cannot continue. Furthermore, he found the singular 
character of the limiting line, such as the infinite acceleration and infimite 
pressure gradient. Th. von Karman (Ref. 7) demonstrated this fact for the 
general two-dimensional flow, He also suggested that the limiting line is 
the envelope of the Mach waves (Fig. 1) and thus can only ocour in supersonic 
region. He also took its appearance as the criterion for breakdown of iser- 
tropic irrotational flow. This renoral two-dimensional theory was established 
later by both iingleb (Ref. 8) and Tollmien (Ref, 9). Tollmien corrected 
some mistakes in Ringleb's paper and in sdditian, showed that the flow de= 
finitely oannot continue beyond the limiting line, The later fact introduced 
a "forbidden region" in the flow bounded by the limiting line. This physical 
absurdity can only be avoided by relaxing the condition of irrotat ionality. 
But as stated previously, for non-visoous fluids, the transition from a flow 
without vorticity to that with vorticity can only be accomplished by shook 


waves, which at the same time also cause an increase in the ontropy. 


However, before one can conclude that the appearance of limiting line, 


or the cnvelopeof Mach waves, is the general condition for breakdown of 
isentropic irrotational flow, one must prove that the singuler behavior of 
limiting lines are general and not limited to twoedimensional flow. This is 
the purpose of the present paper. First the property of limiting line in 
axially syrmetric flow will be inv-stigated in detail. Then the general 
three-dimensional problem will be sketched. These investigations confim 
the results of Ringleb, von Karmn and Tollmien for these more general cases, 

Therefore, by considering only the stesdy flow of nom-visoous fluids, 
the criterion for breakdown of isentropio irrotational flow is the appearance 
of limiting line, However, for the actuel motion of a solid body, the flow 
is neither steady nor noneviscous, Small disturbances always occur and almost 
all reel fluids have appreciable viscosity. The small disturbances in the 
flow introduce the question of stability. In other words, the solution found 
for isentropic irrotational flow my be unstable even before the appearance 
of limiting line, end tends to transform itself to a rotational flow involving 
shock waves at the slightest disturbance. If this is the case, the criteria 
concerns not the limiting line, but the stability limit. This problem has yet 
to be solved, 

The effect of viscosity will be limited to the boundary layer if the 
pressure along the surface in the flow direction never increases too rapidly. 
Then outside the boundary layer the flow is isentropic and irrotatioml, 

If the gradient of pressure is too large, the boundary layer will separate 
from the surface. However, at low velocities su h separation only widens the 


wake of the body and changes the pressure distribution over the body. But 


if the boundary layer separates at a point where the velocity outside the 


“be | $ 


boundary layer is supersonic, a^ditional effects may appear. The flow outside 


the boundary layer in this oase can be regarded approximately as that of a 
solid body not of original contour but of a reg contour inoluding the "dead 
water" rerion orented by the separation, It is then immediately clear that 
the ideal isentronic irrotational flow around this new contour my have a 
limiting line, lienoe, the actual flow then must involve shock waves. In 
other words, the separation of boundary layer in supersonic region may induce 
a shock wave and thus extend its influence far beyond the region of separa= 
tion, Furthermore, the steep adverse pressure gradient across a shock wave 
may accent the separation, This interaction between the separation and the 
shock we ve is frequently observed in experiments. 

The above considerations indicate the possibility of the breakdown of 
isentropic irrotational flow outside the boundary layer even before the 
appearance of limiting line, Therefore, the Mach number of the undisturbed 
flow at which the limiting line appears may be called as the "upper critical 
Mach number." On the other hand, sinee shock wave oan only occur in supere 
sonic flow, the Mach number of the undisturbed flow at which local velocity 
reaches the velocity of sound my be called as the "lower oritionl Mach number", 
The actual critical Mach number for the appearance of shock waves and the 
compressibility burble must lie between these two limits, By oarefully 
designing the contour of the body to avoid the crowding together of Wach 
waves to form an envelope and to eliminate adverse pressure gradient along 


the surface of the body, the compressibility burble can be delayed. 


Axially Symmetrio Flow 
The solution of the exact differential equations for an axially 


syrmetrio isentropic irrotational flow was first given by F. Frankle (Ref. 10) 


The method is developed independently by C. Ferrari (ep, 11). Their method 


applies partioularly to the case of supersonic flow over a body of revolu- 
tion with pointed nose, In this case, the flow at the nose can be approxis 
mated by the well-kuom solution for a cone, From this solution, the differe 
ential equstion is solved step by step using the net of oheruoteristios which 
are real for supersonic velocities, In the following investigation, the 
chief ooneern is not the solution of the partiel differential equation but 
rather the occurrence and the properties of the limiting line in an isentropic 
irrotational flow, The general plan of attack is that of Tollmein (Ref. 9). 
However, hore the caleulation is based on the Legendre transformation of 
velocity potential instead of the stream function. 

If 4 is the magnitude of the velocity, /L the corresponding velocity of 
sound assuming isentrop!o process, f the pressure and x the density of 


fluid, the Bernoulli equation aus 
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$ 


In these equations, the subsoript o denotes quantities corresponding to Sg 


, -É 
a d (3) 


and / is the ratio of specific heats of the fluid. Let the axis of symmetry 
be the x s axis, the distance normal to x - axis be denoted by y, and the 
velocity components along these two direstions be denoted by u and v respeo= 


tively (Fig. 2). The x = y plane is, therefore, a meridian plane. Then the 


kinematical relations of the flow are given by the vorticity equation 


and the continuity equation 

ie a 
Eqs. (1) (2) (3) (4) and (5) together with tho relation 7“: 27” sneoity 
the flow oomletely. 

To simplify the problem, a velocity votential P defined as follows is 
introduced: 
e= Py (6) 
Then Fq. (4) is identieally satisfied and Eq. (5) together with Eqs. (1) 
and (2) gives the equation for f Pr 
wv 
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The characteristics of this differential equation, to be oalled the chr 
acteristios in the physical plane is given by Zi Ela » where A0) 
is determined by the following equation 
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It oan be easily seen from this equation that / is real only when 272. 
Thereforo, the characteristios are real only in supersonic regions of the 
flow. 

The meaning of characteristics in the plysical plane is immenintely 
clear if one oalouintes the relatio between the slope of a characteriatia 
and the slope of a stream line ir the meridian or xy plene, By the definition 
of the function g(x,y), the value of g is zero, or constant, along a ohar- 
EEN Therefore, by writing a quantity evaluated at a certain constant 
value of a parameter with that parameter as a subsoript, the slope of 


characteristic in physical plane is 


_. (9) 
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Along e strenm l/ne, the stream function y% defined by following equations 
is constant: 


WIR a, f- -yfr (10) 


Therefore, the slope of a stream line is 
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(8), (9) e (11) give 
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Send d is the Mach Mis given by B= SEN ES Therefore, Eq. (12) shows 


that the characteristics in physical plane are inclined to the stream lines 
by an angle equal to the "ach angle. Such lines are the wave fronts of in- 
finitesimal disturbaroes and are called Mach waves. In ther words, ohare 
acteristios in physical planes are the Mach waves in that plare, There 
are two families of Mach waves inclined syrmetrioally with respect to each 
stream line. 

If to each pair of values of u and v, there is one pair of values of 
x, y, then x and y oan be considered as funotions of u, v. In other words, 
instead of taking x and y as independent variables, u, v oan be used as in- 
dependent variable. The plane with u, and v as ooordinates is called the 
"nodograph plane." An equation in hodograph plane oorresp nding to Eq. (7) 
can be obtained by means of Legendre's transformation. By writing 


UX vg -y (13) 
it is seen that 


A,94 / Sé WA (14) 


Then Eq. (7) oan be written as 
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The charaat ~istios of "d, (15) are given by f l4,r)-0 were f is the solu- 4 UND 
tion of following difforential equation EE I 
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Eq. (16) shows that the characteristics in hodograph plane depends upon the 


values of the derivatives of A which must be obtained from Eq. (15). In li 

other words, the charactaristios in hodograph plane chan e with the flow 

ard are not a constant set of curves as those in two-dimensional problems, ] | 
To obtain the relation between the characteristios ir physical plane | 

and those in hodogranh plane, it is notioed that Eq. (9) can be rewritten y 


4: (ey mh A (17) 
Then Fq. (8) is ms to 
(- BG + AE l) (8), + li- Etat: ae 
However, in EE Fine wee mes the following SE between the B 


differentinls of x and y and those of u and vt 
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By mens of these relations, "d, (18) can be transformed into an equation for 
(du ) and (dr) a This transformed equation can be simplified by using 


Eq. Se? the final relation is 
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m Therefore, if tho first factor of Eq. (20) is not zero, the variations ke? 
M and (dr, elong a characteristic in physical plane must satisfy the relation 


IR ët ZIL jer = fa E^ f^ (a, Hz) ZE v Z0 
f This is tho sane EN for the variations e Sue along a charactere 
[ isties in hodograph plane as can be seen from Eq. (15) and the folloving rw 
lation obtained fron the definition of 

(dr) : (de d d (22) 
The transformed characteristics of physical plane and the characteristics of 
hodograph plane themselves satisfy then the same first order differential 
equation. Therefore, these two types of curves are the sane. In other words, 
the cherasteristies of hedograph plane are the representation of Nach waves 
in u, v plane. 

The Limit. Line 
Eq. (20) shows that if 
d u 5 sp (28) 

then the transformed differential equation for oharacteristios of physical 
plano, or Mach waves is sat'éfied, Therefore, if thore is a line in the 
hodogr&ph plane, aleng which the values of the derivatives of PA are sudh 
| that Bq. (23) is true, then this line when transferred to the physical plane 
will have its slope equal te that of one family of Mach waves. Such lines are 
called as limiting hodograph in u-v plane and limiting line in physical plane, 
Sinee Mach wavos only ooour in the supersonic regions, it is them evident 
that the limiting line mst appear in these regions, The signifioanee of 
the adjective "limiting" will be made clear as other properties cf sueh lines 


are investigated, 
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Now the 1. estion arises: Cen the limiting hodograph be the characterise 
tio in u-v plane? Along a limiting hodogrenh, Eq. (23) rives 


(4) = = RS Ee = dee Ae Tu des (24) 
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where the subsorint y; denotes the value alo:s a limiting hodograph, Now the 
general differential equation for FE » Eq. (15), is true for the whole u-v 
plane, therefore, the equation is still true by differentiating it with 
respect to u and v, "he results ean be simplified bv using jq. (15) itself 
end Eq. (23). Then at the limiting hodograph, 


[t ge 2) ^d, dr] À tdf = a kol kal E) sj Ar (26a) 


- PUBL ABA, + HIER | : 


vu 
E E MA sa ALS 
[Or &)* X, de ike Se nee «ftt aU" 4, au, 
= (GH) E Apy Aer GP 
Eqs. (24), (25a) and (25b) are the only available equations involving no 


higher derivative than the third, Or the other hand, the slopo of a char» 


aoteristic in hodograph plane can be calculated by Eq. (22), 


CG e + (26) 


This equation together with Eq. (16) gives 


0-2) Tf) d z Aufn pE CL, "fen 


dr. 


A BITTE PIPER Ten. 


Therefore, if the limiting hodograph is a characteristic, then fi E / must 
satisfy Ey. (27). However, a simple caloulation shows that it is not even 
possible to obtain a relation between idi) and other quantities not ins 
volving the third order derivatives of E żġa Henoe, a) does not satisfy 
Eq. (21), In other words, the limiting hodogra h is not a characteristio, 
Transferred to physical plane, this means that the limiting line is not a 
Mach wave, But as shown in previous paragraphs, the limiting line is every- 
where tangent to one femily of Mach waves, Consequently, the limiting line 
must be the onveiope of a family of Mach waves. This property of limiting 


line can be taker. as its physical definitio. 


Liniting locograph end the stream Lines 


At the limiting hodogreph both Eqs. (15) and (23) hold, by eliminating 
/ 
one of the second order derivatives, say Iu » the following relatim is 


obtained 


Av. 2 
d "an use (28) 
E ge ee 
The sign before the redical in iq. (28) can be either positive or negative 
but not both. This relation will be used presently to show that the stream 


lines and one family of characteristics are tangent in uev plane, 


From q. (10), the differential of stream function can be enloulatod 


ERT gn "T 


In this eguetion, y can be replaced by A, according to -qe (14) ax the 


us 


dif^orontinla øy end 4 replacol ly the differential Mend 47 &ooordi:g to 


KH (19). Then 


d= WE [er Te ar tHe Lr) de | (30) 


-15- 


Along a stream line, AN - p s therefore the slope of stream line i: hodograph 
plane is given by 


(££ - E E LP A (31) 


At the » d hodograph, Eq. (23) holds, therefore, Fq. (31) together with 


Eq. (28) rives ot 
Ke 
pz» Aur. | = Ba FT (32) 
duy d "(x = f. L- 


where the i before the radical can be either ER or positive oorrese 
ponding to the sim in Eo. (88). 

On the other hand, the slope of the characteristics in hodograph plane 
is determined by "o, (27). By solving tor JE | and simplifying the result 
with aid of To. (15), 


SG n 2 
de 2 v um e. 2A (33) 
Y EZE L Ar 


The sign before the radical is either positive or negative corresponding to 
the two fanilies of characteristics, By using the positive sign in conjuno- 


tion with vositive sign in Eq, (20), and similarly for the negative sign, 


a) l-4 
MIT E VET 


Eqs. (52) and (34) show that the streem lines and one family of ohiraoterise 


(34) 


tics are tangent to each other at the limiting hodograph. ‘his result is 
the same as that obtained for two-d imerlgi onal flow, (xof. 7, 5, 9), “hese 


equations when compared with de (12) for the elope of “ach waves in physical 


plane yields the irteresting result that the stream lines and one family of 


m | Sinee 


characteristics at limiting hodograph are perpendicular to the oorresponding 


Mach waves at the limiting line. 

il-4 = 
gg V, 

Eq. (32) gives the following equation which holds at the limiting hodograph 


(- Ey «245 Heel m 
This equation can be reduced to more familiar form by introducing the polar 


coordinates in u, v planes 
us 206, v= GE 
where © is the angle between the velocity vestor and É-axis, Then Ta. (36) 


takes the form 
(Wir Bly 70 
This can be roparded as the equivalent to kq. (23) for def ning the limit l 
hodograph. Similar relation exists for two dimensional flow, (ef. 7, 8, 9) 
Along a stream line, the ratio between dr, and (da), is given by 
Eq. (31). Sy substituting this ratio into ^q. (19), the differential ay 


and (e), elo-g a stream line is given as 
idi) = pre li FA 
, ed rm +4 dor 
ee EE 
d -v ky telyr 


(da), (se) 


ty 


At the limiting line, ia. (25) is satisfied, Then Lo, (38) shows that at 


the limiting line, the stream line has a singularity. Or, more plainly, (ty 


=l» 


and (dy at these points are infinitesimals of higher order than (de) 
and (dry « By writing 5 for the distance measured along a stream line, 
Eq. (38) gives immediately 


ee a 
P bf Au Se - cR = 


Similarly, 


mei = —— Vitae an A dur ix 
EO Wa hau e 


Therefore, at the limiting line, the accelerations along a stream line is 
irfinitely large. ‘urthermore, since the pressure grudient Ki A elor.g 
a stream line is 
fly -f44 = pf (a), tv (ty | (a) 
the pressure gradient at the limiting line is also irfinitely large. 
2uoh infinite acceleration and pressure gradieitt lead o:.c to suspeot 
that the fluid is thrown back at the limiting line. In other words, the stream 


lines are doubled back at this line of singularity, To investi ate whether 


2 
this is true, tho chamoter of the relation E AS ie Zoe =O alonga 
Stream line hes to be determined. If the derivative of this expression along 
2 
^ strom line is not zero, then des ra PN has only & simple zero 
at the intersection of the limiting line and the stream line, Tonsemently, 
the "ifferectíiols (dk), and (4 Jp will change sign by nessinr throurh the 
Linitine hodorrenh in uev plane along a stroam lire,  Henoe, the stream lines 
will double baok and forma ousp at the limiting line, The derimtive of 


YA.- NY, along the stream Une oan be nr ted with the 
Ud VV 4r 


— 


Ka E Cf vate a n 
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aid of Eq. (30) 


RER du) E Bd euius 
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The expression on the right of Eq, (42) cannot be reduced to zero by the 


(42) 


available relations which consists of Eq. (23), Eq. (15) and differentiated 
forms of Fq. (15). Therefore, the expression concerned generally only has 
8 simple sro nt tho limiting hodograrh end the strem lires are doubled 
back at the limiting line, It will be shorm later that there is no solution 


nossible beyord the limitinr lire, once, the name limiting line, 


Sinee the limiting line is the envelope of the Mach waves in physical 
plane, it is interesting to see whether there is also envelope for the 
characteristios in hodograph plane. The chareoteristios in u-v plane are 
determined by Lq. (26). The envelope to them can be found by climinating 


GI between iq, (26) and the following equation 


(43) 
er er. lie, 
al Ap “ve [Mu p 
which is obtained by equetirg to zero the Tc derivative of Fo. (26) with 


rospect to AR The result can be simplified by "q. (15), and then it is 


simly 
j= uto? + ae Ge 
a' 4 EL (44) 
is is satisfiod by either 
or 
u^ d ue (46) 


al7e 


The first cardition, Eq. (45), when substituted into Eq. (26) gives 
I) a = E wn 

which shows that the oircle of maximm velocity corresponding toea = 0, is 
the envelope to the characteristics in hodograph plane. The second omdie 
tion, Eq. (46) is the spurious selution, since generally the characteris*io 
at 4 = a does not tangent to the cirole 2 =a. Hence a = O is the only 
envelope. 

The lines of constart velocity in hodograph plane are simoly circles. 
Therefore de) m 
AU 4 De (48) 
By means of this relation and 'q, (19), the slope of the lines of ooistant 


velocity is given as 


/ =H 
EE (49) 
4 7 
/ Phy =A, 
This equation together with Eq. (30) gives the following interesting relation 
/ 
; / 
Gi), =- TEN o 
In other words, a line of constant velocity in physical plane is perpendi- 


cular to the stream line in hodograph plane at corresponding point. 


The Lost Solution 
Throughout the previous calculation, the possibility of using the 
Legendre transformation is assumed, This requires that for each pair of 
values of u, v there is one and only one pair of values of x, y. ‘lowevor, 
it is not always true, it is possible to have a number of points i. Lie 
physical plane havi: the sume value of u and v. If this is the cass, chen 
evidently it is impossible to solve for x and y from the gadr of Twiotions 


u = u (xy), v = v (xy). äathemtiocally, the situation is expressed by 


218. 


saying that the Jaoobin Héi v)/als y) vanishes in the physical plane, Or 
U V, - 4 5 =0 (51) 

However, this is also the condition for a funotional relation between u and 

Vp 9.g., V œan be expressed as a function of u, In other words, u and v 

are not independent, "enee if a solution is "lost" or not included in the 

femilv of solutions allowing Legendre transformation, then for that solution, 
y= v/a) (82) 

It is seen that “q. (51) is then identically satisfied, 


By eliminating (à from the ru Penn we obtain 


ail. 4 ( Het hd Ei 4 /l- LE (53) 
This equation oan be rewritten in the S d form by using Eq. (52) 
A. 
Jie dt) Ee ar jü-E Aa du i ER 


The vortioity equation, Eq. (4) ean be expressed as 
45 y, - 4y=0 (56) 
From Eqs. (54 ) and (55), one can solve for Ü, and 


is 


[t E) 458 E GEI] "7 
f! (1-4) - Se (EJ 4 = ar T 


By differentiating the first of :q. (56) with respect to y, the second with 


respect to x, the following relation is obtained by subtraction! 
Zr U, + ak (67) 


deut y $ 

i = BE (58) 

or NES 
/ 


where Hy) se an undetermined function of y. However, :q. (86) anows that 


Therefore 


for lines of constant values of u where // = 4, (dx) + P. (4g), =? 


së 


GI ee du = constant (89) 


Hence, lines of constant values of u and v are straight lines, This ree 
striction reduces the function fiż) sm Eq. (58) to a numerical constant. Put 
Abg) - K » Eq. (58) is then 


/ = 2 (60) 


Therefore lines of constant value of u and v are radial lires passing through 


the point x = K. Thus the lost solution is nothing but the well-known solution 
for the flow over a conical surface, 

From Eq. (59), it is seen that lines of constent velocity are perpen= 
dicular to the tangent of uev curve at the corresponding points. By substi- 
tuting the value of Z from Eq. (57) into Eq. (56a), a relation between u 
and v is obtained: 

2 A 2 2 2 
rfe UB) Eira BB Eo en 
This is the differential equation for the determining the hodograph represente 
ing the flow over a cone. Fig. 3 shows the hodograph for a cone of 30° semi- 
vertex amÓg o and with a velocity at surface of cone equal to 0.35 o. c is the 
maximum velocity, i.e., the value of j corresponding toa = 0. Fig. 3 is 
drawn from data given by G. I. Taylor and J. W. Macooll (Fef, 12). 

It may well be mentioned here that the lost solution for the axially 
symmetric flow is not limited to supersonic velocity only as is the case for 
twoedimensional flow. In fact, Teylor and Maosoll show that for small for- 
ward velocity of the cone, supersonic velocities cour only just after the 


head shock wave. The velocity decreases as the surface of the cone is 


approached. Finally, it becomes subsonic for points near the surface of the 


cone. Fig. 4 shows a few examples taken from their caloulations (ef, 12). 


The dotted curves in the figure are the Mach waves. The dotted straight 
lines are the boundaries between the supersonic and the subsonio rerions, 
Furthermore, spark photographs of conical shell in actual flight teken by 
Mace 11 (Ref. 13) do not indicate the presence of shook waves in regions of 
flow where such transition from supersonic to subsonic velocities is expected. 
Therefore, at least for this particular type of flow, a smooth transition 


through sonic velocity actually takes place. 


Continuation of Solution Beyond the Limiting Line 
Sinos it is shown ir a previous paragraph that the stream lines are 


generally turned beck at the limiting line, the question arises: Is it 
possible to continue the solution beyond the limiting line? Of course, 
there are two ways of co tinuing the solutions The new solution is joined 
either smoothly to the given solution at the limiting line or with a dis- 
continuity. As shown before, the limiting line is the envelope of one family 
of tho Mach waves, then at every point of this line its direotion differe 
from that of stream line by an angle equal to the Mach angle. ut the "ach 
angle is not zero except at points where the velocity of fluid has reached 
the maximum velocity and the ratio A 20. Therefore, the limiting line 
generally does not coincide with the stream line, an’ the discontinuity at 
the junction of the solution at the limiting line cannot be that of a vortex 
sheet, The only other type of discontinuity is the shock wave, However, 
the angle between the limiting line and flow direction is equal to "ach 
angle. Then &ccor?^inr to the result of the theory of shock wevcs, the dis- 


continuity across such & line vanishes. In other words, there cannot be a 


wäll e 


diseonti. it at the limiting lire, Therefore, it is impossible to joi: a 
new solution at the limiting line with a disomtinuity. 

48 to the second possibility of joining a new solution smoothly at the 
limiting line, it is seen that the flow beyond the limiting line must be 
irrotetional acd isertropic since the limitinr line cannot be a shook wave. 
There are only two types of isentropic irrotational flow, namely, one that 
allows the Legerdre transformation and one that does not, the lost solution. 
Inv stigate the second alternative first. If the solution beyond the limiting 
line belongs to the so-called "lost solution", then since the junction at the 
limiting line must be smooth, the values of u, and v at the limiting line 
must also satisfy the :q. (61). But, the slope GL at limiting line is 
given by Eq. (24). The second derivative Kä will then involve the 
fourth order de, vatives of pi . Besides these expressions, the available 
relations are “qs. (23), (15), (25a) (25b) and three more equations obtained 
by differentiating "on, (25) with respect to u and v. However, it is still 
i-possible for A to satisfy an equation like Eq. (61) where no derivative 
of A appears. Hence, the limiting hodograph does not satisfy the equation 
for lost solution. In other words, the "lost solution” cannot be used to 
continue the flow beyond the limiting line. 

The only remining possibility is to continue the flow smoothly by 
another solution obtainable by Legendre transformation, Smooth continuation 
means that the values of u, v SÉ must be the same at the junction, the 
limiting line. Since shock waves do not appear, isentropio relations still 
hold, The density ¢ is determined by velocity only. The value of u, and 
v are determined by the coordinates in hodograph plane. The »osition of the 
limiting lire i. physical is determined by Ku, y; e Therefore, the 


problem oan be stated as follows:At a oertain given curve 4/// 7/A|in the 


hodograph ;;lane, the limiting hodograph, the values of E. 4, ki are 
given. $ is the parameter along the given curve, It is required to deter 
mine a new solution of the differential equation Eu. (15) with these initial 


valuos, First of all, it is seen that with the given data, the le^t hand 


sides of the following equations are givens 


FUJ- 444 CH 
tt he EE (62b) 


Eet, ZE (63a) 
P = [thy f PW td £n Wey (650) 


By substituting those values into Eq. (15), the second degree terma reduce 


"n d 18 Zil f Pie? dla) [ow Ihr (s4) 
+ [RT JR} 


which is linear in daw - Therefore Fei can be uniquely determined by 


Therefore 


Tq. 15). In other words, with the given data, the second order derivetives 
of | at the given curve alh), vi) can he determined uniquely, in spite 
of the fact that the differential equation (15) is of second degroe. Ke 
Friedrichs and H, Lewy (Ref. 13) have shown that under these circumstances, 
the function | within a region R (Tig. 5) bounded by two charasteristios 
and the given curve is uniquely determined except on additional constant. 
Consequently there ca: be only one solution corresponding b the river datu 
at the limiting hodograph. However, this solution is the very one which 
gives the reverse flow at limiting line. Therefore, it is impossible to 


d 


continue the solution beyond the limiting line even by Legendre transformation. Hr 

Since all three altermtives fail to offer a way of continuing the . 2 
solution, the limiting line is truly an impossible boundary to oross, In E 
other words, the rerion beyond the limiting line is a "forbidden recion", | 


This nhysical absurdity can only be resolved hy the br akdom of isentropic 


irrotational flow some distance ahead of the limiting line. 


General Three Dimensional Flow 


— omms 


The methods used in previous seotions for investigating the axially 


syrmetric flow can be easily extended to the general three dimensional 


TTT N 


— 


oase, In the present section, this investigation will be sketched briefly 


and the results indicated. 


Let the three comoonents of velocity along the three coordinate axes 
x, y and z be denoted by u, v, and w respectively, Then by introducing a 
velocity potential // defined by 


u$ n ege MG (65) 


the differential equation for ff of an isentropic irrotational flow can be 


written as (Hef. 7) l 
2 +y y wtf (66) 
Al Ba t yy the n Letz f, td, o Aen "cf | 
if for every triad of u, v, W, there is only one triad of x, y, Zə then | 
4 
the Legendre transformtion can be used, Thus | 
Y = ux +v tw -9 (67) 
and 
= = = 8 
dou eat dues (68) 


The differential equation for / , Eq. (66), is then transformed into 
a^ [BC- FrrcA-G’rAB-H]= nae FL)» A-G) teo) 
4w (48 - 2) dew (GH-AE) + den (ME lant 


4 
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where the ‘ollowing notations are used 
ġ 2 (70) 
E bho el Fak. Gal, Pte 


By analogy with the axially symetrio case, the limiting hodopraph 
surface is defined as the surface in the u, v, w space, or hodograph space, 


where the following relation holds: 


IA HG 
A= wg. 1 Se (n) 
la Fc 


The prepertios of this limiting hodograph and the corresponding limiting 
surface oan be found by considering the behavior of stream lines and charao- 


teristics at such surfaces, 


From Fq. (68), the differentials of x, y and g oan be written as 


dx = A dus H dy + 6 dw (784) 
Aj = H du + B de o F dw qe) 
Aë (3 éi 4 E det C dw Ge 


Along a stream line, the differentialsa dx, dy and dp must be preportional 


to u, v and w respectively. Thus the equation of a stream line in physical 
ay by _ bey 
Cie ye (13) 


where the subsoript y indicates values taken along the stream line. The 


space is 


equation of a stream line in hodograph spase is obtained by eliminating dx, 


dy and ds from Eq. (73) by Eq. (72). The result is 


— yy, (any CR MN (74) 
ECH hut gu  Autbur dw u+ fu HEW 


where A is the co-factor of A An the determinante A of Eq. (71), Z the 


cofactor Witz. Eq. (74) oan be used in turn to eliminate tw of the 
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three differentials du, dv and dw in the right of Eq. (72). The result is 


Mix) = tia — | 


y Aut dv fu (75a) 
dy) _ r A ar i 
Ay = But Prt dw un 
u — Sy ^ (75e) 
dey, Ju+ {r+ ew 


At the limiting surface, Á =0 as defined by Eq. (71), therefore the stroam 
lines have a singularity there, Similar to the axially symmetric flow, 
the stream lines generally are turned back and form a cusp at this surfaoe, 
The acceleration and the pressure gradient are, of oourse, infinitely large 
at such placos, 

The characteristio surface AX, f, 2) =9 in physical space is deter» 
mined by the equation 
aL fy + by +4, = Wf + A eA drug den fitu f, fy (18) 
Since this equation is a second degree equation, there are two families of 
surfaces passing through each point. These surfaces are the wave fronts of 
infinitesimal distrubances in the flow and can be called the "ach surfaces. 7 
The characteristic surface Aiu, V, ur) = in the hodograph space is deter- 


a tr 
aftBr c) uh HM + URB) MEL LLL HE] cy 
u Ch BL" AE Je nl PO Me] 
«w^ [ Bf, Se Hi 4 Tr Au] Ht ht GAL FATICA 
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By transfo mins Ta. (76) for Mach surfaces to hodograph space, it can be 
shown that tho transformed equation is satisfied by either the ohumoterise 
tios in hodograph space determined by Eq. (77) or the limiting hodograph 
determined by `o. (T1). Therefore, here again the limiting surface is "he 
envelope of a family of “ach surfaces, 

By using “qs. (74) and (77), it is possible to show that the strean 


in n lines ín hodograph space are tangent to the oharaoteristio surfaces at the H 

limiting hodogreph. Furthermore, by using Eqs. (69), (71) and (74), the 2 
inclination of the stream lines at the limiting hodograph oan be calculated. 

In fuot, tf. (45 )~= (du) (du) o [dw]? F= M) i) u the following rela Ki 

WW tion is obtained D 


4 5 
Can = 4 or L (78) — 


E This relation is really equivalent to Eq. (32). In other words, at the limite 
(Ri ing hodograph, the inolinatian of the stream lines and characteristics from 

| the 7 = constant surface is equal to the Mach angle (Fig. 6). It thus sooms 
ja the break down of reneral steady isentropic irrotational flow of non-visoous 
fluid is connected with the appearance of the envelope of Maoh w vos in 

e physical space and the tangency of stream lincs and characteristics in 1 


hodorraph space. 
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Limiting Line as the “nvelope of Mach aves 


Stream Line and Velocity Components in an axially symmetric 
Flow 


Hodorraph of the Flow over a Cone of 300 Helf Vertex Angle 
and a Surface Velocity 4 equal to 0,3BC, 


Flow over Co es of Various Vertex Angle involving Subsonic 
Rerions. (© Half Vertex Angle, 4 = Velooity over the 
Surface of Cone 


Rerion R where the Solution is uniquely determi ed by riven 
Data at the Limiting Hodograph. 


Geometrical Kelations between Stroam Line and Characteristic 
Surf&oe et Limiting Surface in lodorraph 59600, 


7 SUPERSONIC 
REGION 


JUBSONIC 
REGION 


J] -LIMITING LINE 
(ENVELOPE OF WAVES Q) 


REVERSE 
FLOW STREAM 
LINE 


HODOGRAPH PLANE 


— LINE OF CONSTANT U, U 


CHARACTERISTICS 


_ TANGENT TO 
d TREAM LINE 


N 
iv CHARACTERISTIC. SURFACE 
A. THROUGH P de 


S. x 
23, VELOCITY AT P 
\ 


Kee 
& 

2 

| 

Go 

\ 

\ 
= 


SAE 


d 


TA) 


ib 


SON STAN” 


124 Web 
03 04 0.5 
HODOGRAPH PLANE u 


LINE OF CONSTANT U, U 


PHYSICAL PLANE 


TANGENT TO 
STREAM LINE 


P Y-CHARACTERISTIC. SURFACE 
THROUGH P 


— CHARACTERISTICS i $ SR j 
191521 3 VELOCITY AT P 
^ \ 


STREAM LINE 


8-10" "ve-04 
| 


QUARTERLY OF APPLIED MATHEMATICS 


BROWN UNIVERSITY + PROVIDENCE, R. I. 


EN fu y BOARD OF EDITORS: John M. Lessells 
K uU db Hugh L. Dryd MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
^ ugh L. Dryden CAMBRIDGE. MASS. d 
OH NATIONAL gefuer Ivan S. Sokolnikoff A 
i í PAO UNIVERSITY OF WISCONSIN séi n 
8 Le Thornton e. Fry MADISON, Wis. y 


. 
+ BELL TELEPHONE LABORATORIES John L, Synge 
M v 463 West ST., New YORK, N. Y UNIVERSITY OF TORONTO mL. 
nd J TORONTO, CANADA i 


Theodore von Kármán j t 
W. Prager, Mai E 
Sieten, Fee or Een April 12, 1943 ova eg 


ke PASADENA. CALIF. PROVIDENCE. Rt 


Dr. Hsue-Shen Tsien 
California Institute of Technology 
Pasadena, California 


Dear Dr. Tsien: 


I just received the program of the Meeting 
of the Mathematical Society at Stanford University 
on April 24, 1943. I see that you will address this 


meeting on The limiting line in mixed subsonic and Y 
supersonic flows. If you have not yet made other 

now. arrangements concerning the publication of this 
manuscript, I should like to suggest that you consider 


the Quarterly of Applied Mathematics as a possible 
means of publication. 
Very sincerely yours, 


de 


W. Prager i 


WP:G 


UNIVERSITY OF MARYLAND 
CoLizse PARK 


April 24, 1943 


Dr. Hsne-Shen Tsien 
California Institute of Technology 
Pasaderra, California 


Dear Dr. Tsien, 
I have noted tuat you plan to speak upon "The limiting 


line in mixed subsonic and supersonic flows of compressible 
iluids" on April 24th. 


This sounds very interesting to me in the light of 
| some recent investigations I have undertaken, and I hope you A 
will be kind enough to send me a reprint of your lecture when L 
it becomes available. l 


Sincerely yours, v l 
a 


Daan Ti Paa s 


MNES Monroe H. kartin 


Robert R, Dexter 

Institute of “he Aeronautieal Seiences 
1505 RCA Building West 

$0 Rockefeller "ines 

New York, Now York 


Dear Mr, Dexter: 


Fnolosed is my paper entitled "The Limiting Line in Mixed 
Subsonic and Supersonio Flow of Compressible Fluids." I am 47 
TB Submitting it for consideration of publication in the Journal, j 


Sinoerely yours 


H. S. Teien 
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May 13, 1943 


Dr. H. S. Tsien 

Daniel Guggenheim Aeronautical Laboratory 
California Institute of Technology 
Pasadena, Calif, 


Dear Dr. Tsien: 

Thank you for submitting your paper "The 
Limiting Line in Mixed Subsonic and Supersonic Flow 
of Compressible Fluids" for possible publication in 
the "Journal", I am putting the paper through the 
ugual procedure and will notify you of the comments 
of the Editorial Board as soon as possible. 

Best regards 


Sincerely yours, 


Editor 
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, E 
Dr. Hsue-Shen Tsien "ud 
California Institute of Technology M 
Pasadena, California "M 


Dear Dr. Tsien: 


I am writing with regard to your address, THE | 
"LIMITING LINE" IN MIXED SUBSONIC AND SUPERSONIC é 
FLOWS OF COMPRESSIBLE FLUIDS, delivered at the re- * 
cent meeting of the American Mathematical Society L 
at Stanford University. Would you care to offer a 
manuscript for possible publication in the Bulletin? 

I think we could allow about eight to ten printed 
pages. 


As you probably know, papers in applied mathe- 
matics are subject to censorship and, therefore, 
there is some risk of having publication held up 
for an indefinite period of time. I hope, however, \ 
that in spite of this risk I shall receive a favor- $ 
able reply from you. 


Sincerely yours, 


Paul A. Smith 
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July 22, 1943 


Mr. Hsue-shen Tsien 

Guggenheim Laboratory 

California Institute of Technology 
Pasadena, California 


Dear Mr. Tsien: 


Enclosed is your paper, "The 'Limiting Line' in 
Mixed Subsonic and Supersonic Flow of Compressible Fluids", 
which you submitted for publication in the "Journal". This 
paper was returned to us by the War Department with the 
following comment: "Objection is taken to unclassified pub- 
lication of the article entitled, 'The Limiting Line in Mixed 
Subsonic and Supersonic Flow of Compressible Fluids'. This 
paper contains information that would be of aid to the enemy." 


Your paper is an excellent one and was highly 
recommended for early publication by Dr. Dryden. I regret 
that we are unable to publish it at this time. Possibly, if 
you do not find means of publishing it as a classified report, 
you may be willing to submit it at some later date for pub- 
lication in the "Journal". 


Thank you for submitting the paper to us. I hope 
you will continue to send future papers for possible publica- 
tion. I should appreciate your acknowledging the return of 
this material to you for our records. 

Best regards. 


Sincerely yours, 


aussen 
Secretary 
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September 14, 1544. 


Mr, R. G, Robinson 

National Advisory Committee for Aeronautics 
1500 New Hampshire Ave., Dupont Circle 
Washington 25, D. C. 


Dear Mr. Robinson: 


About a year and a half ago, 4 submitted to NACA a paper entitled 
"The Limiting Line in Mixed Subsonic and Supersonic Flow of Compressible 
Fluids" at the suggestion of Dr. von Karman for considerstion of publica- 
tion as a classified report. Six months ago, Dr. von Karman w:8 in 
Washington and has talked about this matter with your staff, und found 
that the paper was mislaid. However, Dr. von Karman told me that it was 
decided to publish this paper at an early date. 


Since another six months have passed, l wonder whether you will per- 
mit me to injuire again about the publication of the paper. âs this paper 
discusses a subject of great current interest, its early public:tion should, 
perhaps, be considered. 


l apologize for troubling you with such a trivial mutter, und 4 remain 


Very sincerely yours, 


Hsue-shen Tsien. 


PU USC n ELS Xu 


JEROME C. MUNSAKER. SC. D., CHAIRMAN RESTRICTED 


LYMAN J. BRIGGS, PH. D., VICE CRAUDMAX 
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REAR ADMIRAL JOHN 5. MCCAIN, U, $. M. 


1500 NEW HAMPSHIRE AVE., DUPONT CIRCLE 


LANGLET MEMORIAL AERORAUTICAL LABORATORT 


AMES AZRORAUTICAL LABORATORY 


LANGLET FIELD, HAMPTON. VA. WASHINGTON 25, D. C. 
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AIRCRAFT ENGINE RESTABCĦ LABORATORT 
CLEVELAND Aurogr, CLEVELAND, O10 


Dr. Hsue-Shen Tsien 
California Institute of Technology 
Pasadena, California 


Dear Dr. Tsien: 


With your letter of August 18, 1943 you forwarded 
to the Committee a copy of your paper "The 'Limiting 
Line' in Mixed Subsonic and Supersonic Flow of Compress- 
ible Fluids" for consideration of its release as a 
classified Committee report. 


Study of this report by the Committee's technical 
staff has been completed and I am pleased to advise you 
that your paper will be published as a Restricted Tech- 
nical Note of the Committee. If it appears advisable 
to refer the report to you in its final form prior to 
publication, and I suspect that this may be the case 
because of the amount of mathematical work included, 
the Committee would like to refer a proof copy to you. 


Thank you for your interest in making the paper 
available to the Committee's staff and in offering it 
for publication. 


Very truly yours, 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


kä 
c—d——-— MONS TRATTE CEST -— d 


G. W. Lewis 
Director of 
Aeronautical Research 


RESTRICTED 
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September 18, 1544, 


Dr. u. Hi, Lewis 

Director of Aeronautical Research 

National Advisory Committee for 4eronautics 
1500 New Hampshire Ave,, Dupont Circle 
Washington 25, D. C, 


Dear Dr, Lewis: 


It is indeed a happy coincidence that only two 
days sgo I wrote to Mr. R. G, Robinson injuiring about the pub- 
lication of my paper "The ‘Limiting Line’ in Mixed Subsonic and 
Supersonic Flow of Compressible Fluids" and then I received your 
letter of September 14, 1944. confirming its publication as a Re- 
stricted Technical Note. Following your suggestion, I would like 
to sec the report in its final form as the paper involved con- 
siderable mathematical work. 


I thank you for informing me, and 1 remain 


Very sincerely yours, 


HsT pl Heue-shen Tsien. 


JEROME C. MUNSAKER, SC. D.. Cmanman RESTKICTED 


LYMAN 1. BRIGGS, PH. D., VICE Cuamsas 
VICE ADMIRAL JOHN $. MCCASIL U. 5. IL 


REAR ADMIRAL E. M. PACE. U. $. L NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


MAL GEM. OLIVER P. ECHOLS, U. S. A. THEODORE P. WRIGHT. SC. D 
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WASHINGTON 25, D. C. 
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Amar? ERGUN RESEARCH LABOBETORT 
CLEVELAND Amos. CLEVELAND, Oo. 


Dr. Hsue-Shen Tsien 
California Institute of Technology 
Pasadena, California 


Dear Dr. Tsien: 


M I now have more complete information on your U d 

idi paper "The 'Limiting Line! in mixed Subsonic and n] 
| Supersonic Flow of Compressible Fluids" and in trans- di 

Wi mitting it I wish to acknowledge your letter of ! 


tal- September 18 and one of September ll, addressed to 
n Mr. Robinson. 


ie 
A mimeographed proof copy will be ready approxi- 
d mately three weeks from now and will be forwarded for 
] your review before proceeding with the duplicating 

l process. The Committee will be pleased to have your 
1 comments at that time and the paper will be published 
1. and released immediately thereafter. 


Very truly yours, 
j 


Ip NATIONAL ADVISORY COMMITTEE 
y FOR AERONAUTI CS 


| Desc | 


1T For G. W. Lewis 
! Director of 
Aeronautical Research 
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Dr. Hsue-Shen Tsien, | 
Celifornia Institute of Technology, 
Pasadena, California, i 


Dear Doctor Tsien: 


With reference to mv letter of September 23, there is forwarded V 
herewith, by registered mail, for final check before release, a proof 1 
copy of the restricted Technical Note No. 961 entitled "The 'Limiting | 
Line' in Mixed Subsonic and Supersonic Flow of Compressible Fluids," 
of which you &re author. 


It is requested that any necessary changes or corrections be 
indicated in red and that desirable changes or corrections be in- 
dicated in black in the enclosed copy of the Technical Note. 


As it is desired to release this Technical Note in the very 
near future, I would appreciate your returning the corrected copy, 
together with the enclosed original material, as promptly as pos- 
sible. 


Very truly yours, 


l Director of 
Aeronautical Research. 


Enc. Proof copy and orig. text and figs. 


RESTRICTED 


October 26, 


Dr. G. W. Lewis 

Director of Aeronautical Research 

Nutional Advisory Committee for Aeronautics 
Washington, D. C. 


Dear Dr. Lewis: 


Follosing your request in the letter of Oc-ober 17, A 
1944, I am returning the corrected proof oup? end origin-l manu- il 
script of my paper entitled "The 'Limiting Line! in Mixed Subsonic Ih 
end Supersonic Flow of Compressible Fluids". There are only two , 
short correctione necessary and they appear on puge 6 and page 24 
of the proof copy. 


It is indeed an honor to publish my paper «s a techni- 
cal note of the Committee and I deeply appreciate it. - 


Very truly yours, 


Hsue-shen Tsien. 


HST/pl 
Enel. 
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